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TI:E production and analysis of MICftOCELLULAR FOAM 

by 

JAN£ ELLEN MARTINI 
Abstract 



A method for producing microcellular thermoplastic foam (cell 
diameter 2 to 10 y) is described. The material produced exhibits tensile 
strengths superior to those of conventional foam. Unfortunately, the 
cycle times are too long for commercial application. To gain insight into 
the design of a more feasible process, the formation and growth of foam 
was analyzed theoretically. The nucleation rate was predicted by class- 
ical nucleation theory applied to a matrix supersaturated with gas. 
Growth of the cells by diffusion ana matrix relaxation was estimated and 
the size of the cells determined as a function of time. 
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Figure l-l Micrograph of Microcellular Foam 



Figure i-i Micrograph of Microcellular Foam 
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I. INTRODUCTION 

As the price of petroleum rises, so does the price of plastic. 
Currently, for high volume plastic p.ru, the material accounts for about 
60% of the total cost (Hardenburg, 1980). In an effort to reduce the 
plastic required fir ;a |>ar|^^ a noy^l fpamini meth^^^^^^^ 
The process vnlT be discussed: in.g,^eate,|aem^^^ thermo- 
plastic sheet is sup^rsat^^^ 
transition t^mp^ra*ui^^^^^^^^^^^^ 

It will ^oam. r^^ ^^^^^ 
2 to 10 u l-v^tameter^arurVxhibi^^^ ^^^^^ 
of conventional foams. Whireas ^thevVtr^^^ ^^^p^ 
to about V5;that of the; pa^^^ 

this foam does not drop more than lis reduction in densi^^^^^^ of 

theparent dpnsity^ this^f^^^ 

more than four times that of conventional foam. 

The microcellular foam shows some other interesting qualities. The 
cells are locally uniform, although those in the center are generally 
larger than the ones near the surface. The surface itself is composed 
of a thin unfoamed skin. This skin probably results from the gas 
diffusing out through the surface rather than forming bubbles. Because 
of the skin, the foam has excellent surface characteristics. In fact, 
by pressing it against a polished plate just after foaming (while it is 
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Still near T^), the surface can be made to shine. Since it Is not porous, 
the surface can be printed and decorated by most of the same techniques 
that are used on the parent material. 

In this work, only Impact polystyrene was foamed and tested quali- 
tatively. However, crystal grade polystyrene and polycarbonate have also 
been foamed by this method using gas. As long as the solubility of the 
gas In the chosen plastic Is high enough, all thermoplastics should be 
amenable to this process. Whether or not the foam produced will be micro- 
cellular depends on many factors. 

Adequate strength and lower cost make microcellular foam attractive 
for many items from calculator housings to tape dispensers. The material 
is well suited for many packaging applications. Its reduced density may 
make It useful for airplane or automobile interiors, and perhaps even 
boat hulls. Measurements have not yet been made of the insulating capa- 
bilities of this foam, but it may also have potential as a novel building 
material. Since this foam is attractive for conmerclal applications, an 
obvious question to ask Is "Can parts be manufactured from it readily?" 
The answer to this is "Not yet." 

Unfortunately, the current methods of manufacturing foam parts are 
not ideal. Cycle times are long, and expansion of the foam Is often 
limited by the mold, which results in large cell sizes. 

Foamed thermoplastic has existed since the 1940's (Benning, 1969). 
Most of it is made either from pellets copolymerized with a volitile gas 
or with chemical or physical blowing agents. Excellent reviews of non- 
proprietary processes are available (Benning, 1969), (Noiseyev, 1960), 
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(Throne, 1979). (Hansen. 1962). (Frisch and Sauders. 1973). There has 
been quite a bit of work done to minimize cell size and make it uniform. 
Generally, this Involves the use of nucleating agents, of which there 
are several categories: 

(a) Fine-sized materials insoluble in the reins. 

(b) Chemical agents that react exothennlcally to produce 
"hot spots". 

(c) Chemical reactants which release gas. 

(d) Additives that partially destroy the cohesive forces 
between the molecules (Benning, 1969). 

Reported cell sizes, however, rarely fall below 0.001 in diameter. One 
exception is an open-celled polyolefin with p diameter cells (Accural, 
produced by Armak, Inc.) reported by Worthy 0978). 

There are three factors particular to foam which profoundly affect 
Us manufacturabllity: 

(1) Foam is not thermally stable. Once an article has been 
foamed, it cannot be raised far above its glass transition 
temperature without affecting foam morphology. Thus, malcing 
major structural changes to the part once it had been foamed 
can pose difficulties. The temperature sensitivity may have 
assets, however, since the material can be returned to its 
unfoamed state by placing it under high temperature when it is 
above Tg. Under these conditions, gas is reabsorbed into 

the polymer. 

(2) When an article is foamed, it expands. It is difficult to 



maintain critical dimensions through a free expansion process. 
However, if foam growth is limited, pressure is exerted on the 
melt which, in turn, causes the gas supersaturation to de- 
crease. As will be shown later, the amount of supersatura- 
tion strongly influences cell size. 
(3) Heat transfer into and out of the foam is slow. If any 

process demands large changes in termperature, the cycle times 
will be long. For the same reason, the surface can be heated 
quickly without changing the temperature of the core. Thus, 
surface modifications may be made without distorting the 
part. 

These characteristics are not peculiar to microcellular foam; they are 
common to foams in general. The manufacturing methods used for conven- 
tionally foamed parts must cope with these factors. To avoid making 
major structural changes in a foamed part, the material is often foamed 
directly into the final shape. Proper dimensions are usually insured 
by limiting foam growth. For example, foamed parts can be injection 
molded by expanding ("cracking") the mold to allow for foam growth. The 
dimensions of the part can be controlled by limiting the mold eApansion. 
Polystyrene foam cups are often made from pentane-filled beads. The beads 
are poured Into a meld and heated. They expand until the mold is 
filled. 

Many of these problems may arise because process designs have 
adhered to the traditional notions of using molds to form parts. The 
niT-Industry Polymer Processing Program is currently developing a 
method for forming parts sequentially, much as sheet metal is formed. 
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This method may be ideally suited for production of parts made from foam, 
since it can accommodate all three requirements. The gross aspects of 
the part may be formed before or as the part is foamed. The part may be 
allowed to foam freely, without regard to dimensional accurany. The 
final diinensions may be brought within tolerance by surface modifications. 

Regretfully, the sequential forming of microcellular parts is not 
yet a reality. The current foaming method still poses a severe handi- 
cap. The gas is introduced into the plastic by diffusion at room temper- 
ature, causing cycle times of about one day. Even raising the temperature 
during diffusion only reduces the cycle time to hours, which is im- 
practical for commercial applications. It appears that gas must be in- 
troduced by a different method. 

To make intelligent guesses about how to modify the process to 
shorten the cycle time, but still produce acceptable foam. It is 
necessary to understand which aspects of the current process are responsi- 
ble for the observed morphology. The remainder of the efforts presented 
here are aimed at that goal. First, tests were made on the effects of 
changing various process parameters on tensile strength. From this data, 
the most favorable morphology could be determined. Then the theory of 
foam nucleation and growth was examined to determine the phencnnena re- 
sponsible for creating such foam. Finally, recommendations are made for 
modifications to the foaming process. Before launching Into these dis- 
cussions, however, the process itself will be presented in greater de- 
tall, and descriptions will be made of other experimental procedures 
used in this investigation. 
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II. EXPERIMENTAL PROCEDURES 
The following procedures were used to make and test mlcrocellular 

foam. 

A. Foaming Process 

Precision extruded, 50 mil thick, impact grade polystyrene 
sheet (Monsanto, high impact, Kodak Roll #684-1, REX #3222) was placed 
in a pressure chamber at room temperature. After the plastic had resided 
in the chamber for the desired length of time, it was removed and heated 
above Tg in an oven. To keep the delay constant, 6 minutes was allowed 
between removal from the pressure chamber and insertion in the oven. 

The pressure chamber was 3-1/2 in (3" id) cylindrical tube of 
stainless steel which was capped at both ends. To allow access to the 
inside of the chamber, one of the caps was sealed with an 0-ring and 
threaded tightening collar. The chamber was fitted with a gas-release 
valve (needle typej, and inlet hose. The inlet hose was steel mesh 
shrouded and rated for 1000 psi. The chamber was also equipped with a 
700 psi burst disk. Gas was supplied to the chamber from a cylinder 
through a high pressure two-stage regulator (Harris #IC3 with Matheson 
gauges P/N 63-3143 and PN 63-3133). Since there was no gas flow once the 
system was sealed, it was assumed that the pressure within the chamber 
was the same as that set on the regulator. Different gases could be 
used by changing the cylinder and adapting the regulator. 

The oven was constructed of two horizontal 1/4" thick aluminum 
plates, 14" wide by 22" long. They were held apart 1/4" by rectangular 
aluminum stock inserted between them along the sides. The oven was 
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closed by an aluminum flange on the back side, but the front was left 
open so the plastic sheet could be inserted. This construction avoided 
the large temperature drops which can occur in conventional ovens when 
the door is opened. Since the oven chamber was only 1/4" high and since 
the other three sides were closed, the temperature was not lowered 
significantly by having the front open. Heat was supplied to the oven 
by a Chromalex barrel heater (PHD lOS. 120 V, 1000 W. R47), which was 
wrapped around the outsides of the plates (see Figure II-2a). To prevent 
edge effects, insulation was installed between the heater and the edge 
of the plates. A West time proportioning temperature controller (Model 
^02) was used with an iron-constantan thermocouple to regulate the 
temperature. With this arrangement, a constant temperature was main- 
tained within 1.5*»F in the 12" X 6" region shown in Figure II-2b. 

For the variable temperature tests which will be discussed later, 
the barrel heater was removed, and strip heaters placed on the top and 
bottom plates at either side. These were set at different temperatures 
and controlled separately. A series of 5 thermocouples installed be- 
tween the two heaters showed that the temperature varied linearly with 
distance. The oven plates were uninsulated for the constant tempera- 
ture tests (to prevent damage to the barrel heater) and insulated with 
fiberglass wool for the variable temperature tests. To facilitate 
moving the plastic in and out of the oven, the samples were placed 
between two ..02 in highly polished aluminum sheets. These could slide 
easily in and out of the oven. Mold release (4 Release Agent. Contour 
Chemical Co., Woburn, Mass.) was sprayed on the sheets to prevent the 
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plastic from sticking and buckling as it expanded. 

B. Variation of Process Parameters 

To gain insight into the phenomena involved in the formation of 
foam observations were made of the effects of changing various process 
parameters on the morphology and strength of the foam produced. 

a) Foaming at Differing Temperatures 

Ten samples 3.8 cm X 30.5 cm (1-1/2" X 12") were put together in 
the pressure chamber with (350 psi) Ng gas for 24 hours. Then, after 6 
minutes delay, they were placed side by side in the variable temperature 
oven described earlier for 7.5 minutes. The oven temperature varied 
linearly from 88«C (190«'F) to 116«C (240»C). The foaming temperature 
was considered to be the average across any given sample. 

b) Foaming for Varying Lengths of Time 

Ten 7.6 cm X 15.2 cm (3" X 6") samples were placed In the pressure 
chamber with gas at 350 psi. They were removed after 24 hours and, 
after 6 minutes delay, placed in the oven at 104"C (220-F). Two 
samples were removed after each: 3 min; 5 min; 7.5 mi n; 10 mins and 
15 mins. 

c) Varying Durations of Pressurizatlon and Gas Type 

Two 7.6 cm X 15.2 cm (3" X 6") samples were placed in the pressure 
chamber with gas for 3 hrs at (350 psi). They were removed and 
after a 6 min delay, put in the oven at 105*C v220*F) for 7.5 mins. The 
process was repeated for various chamber residence times up to 116 hrs 
and with N^, COp and the gas. 
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(a) 




putiine of barrel heater 



A 



•-constant femperatura area 
Top View 
(b) 

Figure INI Constant Temperature Oven 



Figure iN2 Pressure Chamber 



hiGURE ii-2 Pressure Chamber 
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Figure I 1-3 Measurement of Cell Diameter 



C. Characterization of Foam Morphology 

The foams were characterized by their cell size distribution and 
density. 

a) Cell Size Distribution 

A scanning electron microscopy was used to determine cell size 
distribution. The foam was fractured by bending it sharply. It was 
mounted on end in a holder (beeswax seemed to be a successful mounting 
medium), and sputtered with gold to form a conductive surface. The 
sample was then observed with an cil«clron fn*croac/»^. 
Micrographs were made of regions of the mid-thickness of the sample. The 
cell size distributions were found by measuring the diameter of each cell 
in the lower right quadrant of the micrograph, then plotting the number 
of cells In each size range. If the cells were not sperical, the measure 
ment was made 45" to the major and minor axles, as shown In Figure II-3. 
It was assumed that the material fractured through the center of the cell: 

b) Density Measurements 

Since some of the foams were quite light. It was difficult to 
use the available density measurement system based on fluids. However, 
the tensile test specimens were all machined from the same template, so 
relative density could be found with reasonable reliability by weighing 
the test specimens and dividing by their thickness. Call this quantity 
£. Changes In thickness across the sample were generally negligible and 
even in extreme cases less than 2%. The actual density can be found by 
multiplying the density of the parent material by the ratio of g_ foam 
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to p unfoamed. 



p foamed = p unfoamed 



£ foam 



£ unfoamed 



D. Determination of Tensile Strength 

Tensile tests were made in accordance with ASTM specification D 638 
(1979, Part 35). The tests were done on Type IV (4.5 in long) specimens. 
See Figure II. 4 for dimensions. They were pulled on an Instron (model 
1125) fitted with a D-cell at a rate of .2in/min at room temperature. 
The thickness of the samples were measured with a 1 in micrometer before 
they were pulled. 
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III. STRENGTH OF FOAMS 



To determine which factors contribute to the strength of micro- 
cellular foams, tensile strengths have been measured for foams produced 
under a variety of processing conditions. Due to limitations in the 
amount of material available, each experiment was done once, and was not 
repeated. The results are shown here to give an indication of the trends 
exhibited. They should be repeated several times to estimate and ira- 
prove accuracy. 

^' The Effect of Changing Gas Concentratinn 

To attain variation in the gas concentration in the polymer, the 
material was left in the pressure chamber for various lengths of time 
(see Chapter II for a description of the process and the experimental 
procedures). Cell size is plotted against concentration (or time in the 
pressure chamber) in Figure II I-l. and the strength/density ratio is plotted 
against time in the pressure chamber in Figure 111-2. Figure III-3 
suggests that larger celled foams are weaker. The primary recommendation 
for further improvement in foam strength, therefore, is to reduce the 
cell size even more. As concentration increases, density decreases 
(Figure III-4). The strength/density ratio Is affected by a trade-off 
of density with cell size. Although a trend between the ratio and con- 
centration is not apparent from the data shown in Figure III-5 it should 
be noted that the ratio is greater than that for the parent material. 
Here, especially, the experiments should be repeated to clarify this 
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Figure III-l Cell Size Vs, Gas Concentration 



30 - 




Xt!8uap/(||6u84S 



i 

Figure I n-2 Strength/Density Vs. Pressure Chamber Time 
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Figure III-^i Density Vs, Gas Concentration 
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Figure III-5 Strength/Density Vs. Cell Size 
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data since It will point to an ideal concentration, 

B« The Effect of Changing Foaming Temperature 

As the temperature at which the foaming takes place increases, the 

density decreases (Figure III-6), the cell size increases (Figure III-7) 

and the strength decreases (Figure III-8). Discussion of the reasons for 
this are found in Chapter V. Section A. 

C. The Effect of Changing Foaming Time 

The data, shown in Figure I I 1-9, has so much scatter that trends of 
cell size versus oven time are not distinguishable. From work described 
later (Chapter IV, Section C), the cells should be expected to grow about 
3| over the period from 3 to 15 minutes (calculated from the graph 
IV-21). The strength of the plastic, however, does decrease with the 
length of time the plastic is left in the oven ^see Figures III-IO and 
111*11). It Is believed that this reduction in strength reflects the 
loss of orientation imposed around the cell by nucleation and diffusional 
growth. 
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Figure III-6 Density Vs, Foaming Oven Temp 
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Figure III-7 Cell Size Vs. Foaming Oven Temp 
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Figure III-8 Strength Vs. Foaming Oven Temp 
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Figure II 1-9 Cell Size Vs. Time in Foaming Oven 
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Figure III-IO Strength/Density Vs. Time in Foaming Oven 
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Figure III-U Density Vs. Time in Foaming Oven 
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IV. THERMODYNAMICS OF FOAM FORMATION 

Introduction 

The formation of a foam cell occurs in two stages: nucleatlon and 
growth. During the nucleatlon stage, the bubble embryo grows from a region 
of high gas concentration to a stable cell having distinct boundaries. The 
cell then grows by diffusion of the gas from the surrounding supersaturated 
matrix, and finally by relaxation of the matrix. 

To enable a cell to form, the decrease of supersaturation must supply 
the energy needed to form the cell and effect the accompanying entropy rise. 
Obeying the second law of thermodynamics, the entropy of an Isolated system 
at constant temperature must increase above dQ/T during an irreversible 
process. The excess energy contributed to entropy is a criterion for de- 
tennlnlng the spontaneity of the process; the greater the energy excess, 
the more readily the reaction can take place. Classical nucleatlon theory 
Is based on this concept. Under highly irreversible conditions, when large 
amounts of excess energy are usurped by entropy, bubbles form quite rapidly. 
If the process is reversible, bubbles will be nucleated at a rate of zero 
(i.e., they will collapse as quickly as they are formed). If not enough 
energy is available to produce the necessary entropy rise, bubbles will not 
form unless energy is added to the system. 

A. Determination of Gibb's Free Energ y 

The amount of excess energy contributed to entropy can be calculated 
by combining first and second laws of thermodynamics. 
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TAS - A£ - AW = excess en tropic energy >^ 0 



(IV-1) 



For a system at constant temperature and pressure, and for which there is 
no shaft work (so AW = PaV. pressure times the change in volume), the 
excess entropic energy just described is the negative of the change in 
61bb*s free energy (G). 



AE + PAV - TAi = A6)t d 

'o* % 



(IV-2) 



G is defined by 



G = U + PV - TS 



(IV-3) 



Thus, for a reaction to occur without additional energy input, AG must be 
negative. For it to occur rapidly, AG must be large and negative. 

The conditions on which AG depends are the concentration of the 
gas (X), tjje solubility of the gas in the matrix (K), the surface energy 
coeffioleni (y), the dynamic modulus (E) of the matrix material, tempera- 
ture (T) and the radius of the bubble (R) after nucleation is complete. 
Dimfensional analysis yields units of energy from the following combinations 
of those parameters: 



KTRyR 
KTRER^ 



T • RR^ 
2 



ER* 
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Where R is the universal gas coefficient. Thus, for large bubble ridii, 
contributions from concentration, solubility and dynamic modulus are ex- 
pected to dominate the expression for AG. For small bubble radii, the 
surface energy coefficient and solubility will be more Important. 

AG will be calculated to understand the contributions of each of 
these factors more explicitly. The system to be used for this calculation 
is a volume of solution (V^) and the gas it contains. The system is as- 
sumed both adiabatic and Isolated from gas transfer. The volume, V^, Is 
such that it supplies all of the gas that bubble needs, but is not larger 
than necessary to do this. It is assumed that the system will change 
quasi statically from an unfoamed state (state 1) to a state containing one 
gas bubble with distinct surface boundaries (state 2).. The following 
assumptions will also be made: 

In state 1: 

A) The solution Is uniform and ideal (i.e., it follows Henry's 
law). 

B) The temperature and pressure are uniform throughout the 
system. 

In state 2: 

A) Even though gas has left the solution matrix, the matrix 
is Incompressible, so the volume of state 2 Is plus V^^, 
where V^^ is the volume of the bubble. 

B) The gas contained In the bubble Is Ideal. The matrix sur- 
rounding the bubble is still a uniform, ideal solution (al- 
beit with a lower concentration of gas). 

C) The syst^ is adiabatic and temperature change is negligible 
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(temperature is therefore still uniform). 
D) Pressure Is uniform in the solution surrounding the bubble. 
One final assumption bears some discussion. It was assumed that the 
nucelation process takes place so rapidly that no gas diffuses to the 
nucleation site. Such an assumption means that nucleation consists only 
of forcing the matrix material from the region which will ultimately be the 
bubble, leaving behind the gas which will fill the bubble. Thus, the 
volume of polymer forced from the bubble region had to contain enough 
gas to fill the bubble (plus any residual gas which was ushered out with 
the polymer). This means that the initial volume of the system, chosen to 
be only that plastic which supplies gas to the bubble, is that of the 
final bubble size. The final volume. V^. of the system is therefore twice 
that of the bubble size, or 8/3 ttrI Note that the assumption that no 
diffusion takes place means that no concentration gradient develops in the 
exiting plastic. 

The minimum concentration necessary for bubbl3 nucleation is 

X - Nb . y 

^1 Vr*^R (IV-4) 



Where is the number of moles of gas molecules in the bubble, and X„ i 
the residual gas concentration. Since the gas was assumed ideal. 



R -5 



M a _b b 

"b rY (IV-5) 



By Henry's law. « P^^K. where K is a constant detennined in Appendix E. 
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(IV-6) 



Henry's law has been experimentally shown to be a good approximation for 
the polymer-gas systems considered in this work (Frisch, 1965). 

The change in Gibb's free energy for the process is the difference 
between the free energy of the final and initial states. Calculating the 
Gibb's free energy for the first state is fairly straightforward. The 
energy stored in the system, E^, is simply the internal energy, U. For an 
ideal solution of two components, 

2 

ii = U, = I y,.n. - PV + TS (IV-7) 
' ' i = 1 ^ ^ 

where \i. is the chemical potential of the ith component and n^ is the 

number of molecules of that component (Denbigh, 1968; Katsopoulos and 

Keenan, 1965). Since the PV and TS terms cancel when this is substituted 

into the definition of Gibb's free energy, 
2 

= 2 M^n. (IV-8) 

For the second state, the Gibb's free energy is most easily calcu-r 
lated by realizing that it is the sum of two parts, the free energy of the 
solution (including the surface of the bubble), 6^, and the free energy of 
the gas. 6g. 

The energy, E^, stored in the solution is the internal energy of 
the solution plus the energy stored in the elastic deformation (D) plus 
the surface energy (S). So 
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2 

§5 ^ Viini - PV + TS + D + S (IV-9) 
i = 1 

Thus, the Gibb's free energy accommodated in the solution is 



S = W*^^^^ (IV-10) 



Where the subscript 2 denotes state 2. 

The energy contained in the gas in the bubble is the same as that of 
a mixture with one component. 

E = Ugnj^- PV + TS (IV-n) 

where n^^ is the number of gas molecules in the bubble, and y is the 
chemical potential of the gas. For an ideal gas, the chemical potential 
is the potential at a reference pressure (p^) plus the potential arising 
from the difference between the actual pressure and the reference pressure 
(RTlnp), 



Pg = P° + RTlnp (IV-12) 



Thus, combining equations 11 and 12, the Gibb's free energy for the gas ^n 
the bubble Is, 



^ v\ + nbRTlnp^ 
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The total Gibb's free energy for the system at the second state is 
therefore. 



The change in 6 for the process is the difference between G., and 6«: 



Since the solutions follow Henry's law, it is possible to substitute 
more meaningful expressions for the chemical potential, namely, express 
U in terms of concentration. 

First, the summation terms can be simplified by making the assump- 
tion that the gas concentration is low. Thus, as the gas leaves, the con- 
centration of polystyrene does not change significantly. Since the number 
of polystyrene molecules remain* constant, hpjWpj Is the same for state one 
and state two. Those terms ca'^cel, leaving 




(IV-13) 




2 

I 

i « 1 




(IV-14) 



' y n. + n„ + D + S - y« n^ + n. RTlnp. 
b gg ^2 9l 9i ° ^ 



(lV-15) 



For any component of an Ideal solution, the chemical potential Is tht 

potential of that component alone plus the potential gained du« AojMliil 
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M = y* + RTlnX 

where X is the concentration of the gas, and y* is the chemical potential 
of the pure component. So, 

AG = + n^jRTlnpjj + ng^{p* + RTlnXg) - v^^ y* + RTlnX^ + D + S 

(IV-16) 

This can be mai:3 simpler still. Since the solution is ideal, there 
is no interaction between the gas and the plastic. The gas will behave 
as though it were in a mixture of ideal gases, and its chemical potential 
can be expressed as that of one component of an ideal gas mixture. 

y = + RTlnp. 

(where p^ Is the partial pressure of the gas). Therefore, A6 can also 
be expressed as 

AG - y\ + njjRTlnpjj + n^^{]i° + RTlnp^) - n^^Cy® + RTlnp^) 

+ D + S (IV-17) 



I. Note that KTlnX is negative since X < l, so the chemical potential 
actually decreases when substances are mixed. This is because entropy 
must Increase. 
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Where is the partial pressure of the gas at state 1. 
Since = Hg^+n,, 

AG = "g^Rtdnpjj - Inp^) + 0 + S 

^ ng^RTln(pj^/p^) + D + S (IV.18) 
Henry's law equates p with KX, so 

AG = ng^.^^TlnCXg/X^) + D + S (IV.19) 

Only those conditions which force this expression to be negative will 
allow nucleation to occur. Many of the quantities in the above expression 
depend upon the final radius of the nucleated cell. Are there certain 
radii which will not be nucleated? 

In order to answer this question, it is necessary to express the de- 
pendent variables in terms of R, All of the quantities except T. K, R. 
and X. depend upon R. T is assumed in variant, K is a material property 
R is the universal gas constant, and X,is an initial condition. 

The dependency of n^. ng^. ng. and can be determined using three 
conditions: 

1. The number of gas molecules in the system remains fixed. 

2. Henry's law applies to the system. 

3. The work done by the bubble to the surrounding solution is 
/PjjdV. 
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From condition 1, the number of gas molecules in bubble, n^^ is the 

difference between the number of molecules residing in the final and 

Initial solutions, n and n . These quantities. n„ and n are re- 

^2 ^1 92 9-i 

lated to the concentrations, and Xg by. 



^1 ' \^^o (IV-20a) 
^2 ' \^^o (IV-20b) 

where is the volume of the polymer solution which supplies gas to the 
bubble (an assumption is made that the concentration Is uniform through- 
out this volume). Since the solution was assumed ideal, 

h ' (IV-21) 

Where p^ Is the pressure inside the bubble. Thus, expressions for 
%* "g^» ^2 °^ ^ follows: 

^2 ° (IV-22a) 
"g, ' \ ' ^1 (IV-22b) 
"92 ^ ^0 * ^^^^ ^IV-22c) 



"b " W^-^i) 



(IV-22d) 
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To calculate Pj^, consider the work done by the gas in moving the 
boundary of the bubble. This work, /p^dV. goes into deformation plus 
surface energy: 



[ Pb^V = D+S (IV.23) 

Where V^^ is the volume of the bubble. /p^^dV can be expressed in terms of 
R Since V = j vR^. 



'O 



AirR'-dR = D + S (IV-24) 



Now, differentiating with respect to R^, and realizing that S = 4TrR^Y, 
^=8irRy 



n . 1 dD . 2y , . 

Pb - ^ HR * ^ (IV-25) 

To calculate the change of the energy of deformation with R, ^- , 
use the standard formula. 



dD = 



3 3 

^ 1 i = 1 ^ij^Sj^V (IV-26) 



Continuity and the assumed incompressibility of the plastic will be uspd 
to determine strain as a function of R. The plastic will be modeled as 
a Maxwell -.olid to determine stresses from the expressions for strain. 
These expressions will be substituted into the above formula so it can 
be evaluated. 



. The plastic is above Tg^ so it should be modeled as a visco-elasti 
plastic solid. The mechanical analog of one example of such a solid is 
depicted 




"^fricrionetement 



However, for simplicity, the material will be modeled as a Maxwell 
solid, . 



In this model, at the instant of loading, the energy is stored elasticall 
then is dissipated through viscous deformation. Thus, the work done can 
be approximated by just considering the elastic part of the model (the 
spring). Although the pressure within the bubble will decrease as the 
strained plastic surrounding the bubble relaxes, the approximation will 
be reasonable if the inflation takes place rapidly. So, the str..ses wil 
be approximated by those of an elastic solid using a dynamic modulus. E, 
which is assumed to depend only upon temperature. So. 



^11 



: J^U ' T~ • (IV-27) 

^ The generalized expressions for strain in spherical coordinates 
are 
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^ s —J. 

'rr ar 



'66 



're 



^inS ^8 



1 ^"a 
* r 3^ ^ r 

+ 7^ + ;r- cot* 
r r 



- 1 f-_L_ !i!r ^ !::e vi 

' 2 Ir,, , 38 3r - r j 



sin(t> 



(IV.28a) 



(IV-28b) 



(IV-28C) 



(IV-28d) 



(IV-28e) 



As the bubble expands, however, the only movement is in the radial direc- 
tion, so = Ug = 0, and the above expressions reduce to 



rr 



3r 



r ^8 



(IV-29a, b) 



(IV-29c) 



Substituting these expressions back into that for dD, 



dD = 



^"r 

W rr^" 



~ Ede««dV 
r 68 



Ede..dV {IV-30) 

J , WW J , ^ 

V \ \ 

The strains can be evaluated in terms of R by invoking continuity 
and the assumption that the plastic is incompressible. The volume of 
plastic which was contained in the bubble must move outward. That volume 
must cross a control sphere drawn at any radius, r. From the geometry 



shown In Figure IV-1 , the volume be expressed as 
r + u. 



^ r 



(IV-31a) 



or 



3 3- 



(IV-31b) 



(Note that R Is the final bubble radius whereas r is the radial coordina 
a variable.) The resulting equation. 



0 ■ 3u^r^ + 3u? r + (u^ R"*) 
r r r 



can be solved for u^ by using the cubic formula. After the ensuing 
algebra, and discarding the imaginary terms. 



= (r^ + R^)^/^ - 



From this result, the strains are 



rr 



R3.-2/3 



- 1 



(IV-32a) 



1/3 



- 1 



(IV-32b) 



de^^ = - 2 
rr 



{I7-32C) 
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Figure IV-1 Geometry for Calculating Strain Imposed 

By Spherical Growth 



de^. - de 
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dr 



(IV-32d) 



These can be substituted Into the expression for dD. 





* 










f 3 1 


-5/3 J, 

r 


dD = 






^*^ -1 

r * 


E 


^2 


r^ 



dR 



(IV-33) 



+ 2 



r p3>l/3 



-1 



1 - 



-2/3 2 
„3 



dR 



Integrating by numerical methods (the program Is shown In Appendix F), 



dD = 12R^Edr 



(IV-34) 



The pressure inside the bubble can therefore be evaluated as 



^ _ 12ER^ . 2y _ 3E . 2y 



(IV-35) 



Substituting this to obtain n^^, and n^. 
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(IV-36a) 



"92= ^ofil^li^] (IV-36b) 



Recall that V^^ was defined as the volume of the polymer solution which 
supplies gas to the bubble. Since it was assumed that no gas diffuses to 
the bubble embryo during nucleation, this volume is | wR^. So, 

"g^ = 1-"^% (IV-37a) 



= 4 erVk + I^yR^ (IV-37b) 



= 4 Er3/K + I - I irR^X, (IV.37c) 



The two other quantities which depend on R, the surface energy, S, 
and the deformatlonal energy, 0, have been evaluated during the course 
of this discussion. 
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S « 4ttR^y 



(IV-38) 




= 4eR^ 



(IV-39) 



All of these values can now be substituted into the Expression for AG to 
obtain an expression which depends on R. 



The values, X, K, y and E» which do not depend on R are ill controlled by 
the foaming conditions. To determine under what conditions foam can be 
produced, AG was evaluated for incremental values of R and various values 



systems at MIT. (The program is called PROG and is in Appendix D.) 

Typical AG versus R curves are shown in Figure IV-3. For 2 moles/c 
AG becomes negative at 1.45X 10"^ y (5.7 X 10"^ in). In order to 
satisfy the energy balance without obtaining energy from external sources, 
the cells must nucleate to this raduis. However, AG goes through a 
metastable equilibrium point at 1.05X 10"^ p (4.13x10'^ in). If 
a cluster of gas i^olecules should occur with a raduis greater than this, 
say 1. iX 10'^u{even if it had to gain energy from the environ.nent to 
reach this configuratioi^, they will be stable and continue to grow. 



AG = ^irR^X, RT In 




+ f ER^ + 4tiR^y 



(IV-40) 



of Xp Y and E. The calculation was made in APL on the MULTICS 




Figure IV-2 Change in Gibb's Free Energy Vs. Radius 

OF Bubble Hmbryo 
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"■3 

Thus, 2.16 X TO M is the minimuin cell size attainable from these con- 
ditions; iall snialler embryos will decrease in size until they are com- 
pletely reabsorbed. 

Wh^t happens if the conditions change? In Figure IV-3, the AG 
versus R curves have been plotted for increasing gas concentrations. 
Since the energy content of the gas per mole increases as the concentra- 
tion increases, the amount of gas needed to overcome the surface ten- 
sion decreases. Consequently, the critical radius decreases with in- 
creasing concentration. This is plotted in Figure !V-4. Since the 
energy released by the gas is a voluiaetric quantity, roughly propor- 
tional to R (this is only approximate since the final concentration is 
a function of R), and the surface energy Is proportional to R^, the 
curve is a hyperbola. This is shown more clearly on the log- log plot 
in Figure IV-5. 

Also involved in ;he energy balance (and hence the expression for 
A6) is the deformational energy. This is also a volumetric term, as 
shown by the dimensional analysis discussed earlier. This term sub- 
tracts from the energy gained by reducing the gas concentration, and the 

resulting quantity must overcorue the increase in surface energy. If 
3 3 2 

CR , OR and SR are the energies associated with concentration charge, 
deformation, and surface formation, respectively, then 

(C - - SR^ > 0 (IV-41) 

In the plot of critical radius versus concentration, only C was allowed 




0 0.5 1.0 \3 

Radius of bubble embryo (A) 



IGURE IV-3 Change in Gibb's Free Energy Vs. Bubble Embryo Radius 

For Varying Gas Concentrations 




Figure IV-^ Critical Radius Vs. Concentration 
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to change. If C is held constant, instead, and D is increased, then the 
energy gained per mole of gas evolved decreases. Thus, the critical 
radius must grow larger and larger until finally, D equals C, and foam 
can no longer be produced. A plot of this is in Figure IV-6. This is 
why cells do not form until the supersaturated plastic is heated. 

B. Production of Cell Size Distributions 

During the proceeding discussion, it was shown that nucleated 
cells must be larger than some minimum size. This size was determined 
for various conditions. Of more interest however, is the size to which 
these cells grow after they have been nucleated and the distribution the 
sizes will follow. 

Bubbles may be nucleated so long as the appropriate conditions 
prevail. If the conditions prevail for five minutes, then the bubbles 
formed initially will have five minutes more diffusion than those formed 
at the end of the time. To determine the cell size distribution, there- 
for^ it is necessary to determine the rate of nucleation and how much 
gas diffuses to the cells per unit time (their growth rate). 

a) The Rate of Nucleation 

In order to determine the nucleation rate, one must rely on 
statistical mechanics. Homogeneous nucleation theory was first estab- 
lished for condensation of droplets from the vapor phase, and then ex- 
tended to other systems. The original theory was developed by Volmer 
(1939), Farkas (1922). Becker and Doring (1935) and Zeldovich (1943). 
They based their work on Gibb's realization that the spontaneity of a 
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process is proportional to exp (-AG/kT). Excellent reviews of the 
development of nucleation theory have been made by Zettlemoyer (1969), 
Blander and Katz (1975), Cole (1974), Reid (1978), and Russell (1980). 

The crucial question in determining the rate of nucleation is 
"how Often does a cluster of gas molecules grow large enough to exceed 
the critical radius?" Following the derivation of Zeldovich (1949) to 
determine this, one must look at the rate nuclei! grow as a function 
of size. The net rate at which bubble embryos one molecule smaller than 
the critical size gain that last molecule is the nucleation rate. The 
embryo growth rate can be expressed as the difference between the rate 
molecules join the cluster and the rate they leave. If the joining 
rate per unit surface is 3, and the leaving rate y, f is the distribution 
of nuclei sizes^and n is the size of a critical cell, the rate of cell 
nucleation, J, equals 

J = f(n-l) 3(n-l) S(n-l) - f(n) y(t)) S(n) (iv-A2) 

where S(n) is the surface area of a cell of size n. 

For a dilute substance dissolved in a liquid, Reiss (1950) has 
shown that 

where V is the diffusion rate and is the mean jump length, and can 
be approximated by the distance between equilibrium positions in the 
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liquid lattice. In the case of a polymer, this would be the mean dis- 
tance between adjacent molecules (see Figure IV-7). Assuming that the 

^tyrene units are uniformly distributed throughout the volume, each 

-22 

polystyrene unit has 1.23 X 10 ml associated with it. If they are 
uniformly packed, each unit is in the center of a sphere, and the dis- 
tance between units is 6.9 X 10 cm. This will be slightly lower 
than the distance between polymer chains since the units will be 
packed more tightly along the molecular axis. However, the bulkiness 
of the unit was not taken into consideration. Neither B S , N nor 

CCD 

Z need to be known with great accuracy since the exponential factor 
will dominate the equation. 3 is therefore approximately 



V is the diffusion coefficient and is discussed 1n Appendix F. X Is 
the concentration of the gas in the polymer and is a variable. 

The rate at which gas molecules leave the embryo can be determined 
for equilibrium conditions when the nucleus does not grow. It is as- 
sumed that this rate remains the same independent of whether an equili- 
brium number of bubbles exists or not. Still following the derivation 
of Zeldovich (1949), at equilbrium. 




6.91 X 10"** cm 



V 



. X 
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J - 0 =. N(ti-1) 3(n-l) S(n-l) 



N(n) Y(n) S(n) 



(IV-44) 



where N(n) Is the equilibrium number of nucleil of size n. Thus, 
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Y(n) = 



N(n-1) S(n-1) eCn-n 
N(n) S(T)) 



(IV-45) 



Using this value for non-equilibrium situations. 



0 = F(n-l) g(D-l) S(n-1) - f(T)) '^^"'"^^(^j'^il^)^^'^'^^ S(n) (IV-45a) 



= N(t)-1) P(ti-1) S(n-l) 



f(n-l) ffn) 



(IV-46b) 



If r » 1, this may be regarded as a derivitive. 



i_ £ 

dn In 



(IV-47) 



where the subscript c indicates that the quantity be that at the critical 
radius. 

As equilibrium, the number of nuclii of critical size is given 

by 



N(n) = 



% exp 



'-A6(t))' 



KT 



XlV-48) 



where is the number density of the available nucleation sites in the 
liquid. 

Writing in a more explicit form, one obtains a differential 
equation 
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Nq, K, T, D and \5.\ are all constants. /iG, and depend upon 

X, and were calculated already during the determination of the critical 
radius. 

Equation 47 can be solved by integration over all embryo sizes 
from 0 to infinity. At very small sizes, f is very close to equilibrium 
so f/N =1. At large sizes, f is small compared to N, so f/N = 0. The 
integration of Equation IV-47 is therefore 



0 1 



(IV-50), 
(IV-51) 



Assuming that the rate of nucTeation is independent from the similarity 
of the actual cell distribution to the equilibrium distribution, 

^ can be solved numerically. More insight will be gained, 
however, by approximation using a Taylor's expansion for AG. Since 
the integrand, l/BSN^ exp l^rpj , has a narrow peak centered around 
(i.e., around R), an expansion about r\ is satisfactory. 

,2 



AG = 



(IV-52) 



(The l^^j term is zero since AG is maximum at the critical point). 
Likewise, 0 and S can be approximated by their critical values. Moving 
what is constant in front of the integral sign. 
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dn . (IV-53) 
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Because AG is a maximum at the critical point. 



t 2 ^ 



3n 



IS nega- 



tive. The integral therefore involves a gamma function: 



2 2 
e t dt 



_ 1 



2a 



(IY-54) 



where 
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r 2 1 
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and t = n - t) . 

So the nucleation rate is 
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J = 2e^S^N„ exp 
ceo ^ 
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As stated before, 3 is the rate gas molecules strike the bubble ■ 
surface, 



ft - 3 D y _ 3 D 

P - y A - ■K- ' g 

^ \i\ ^ 6.91 X 10"" cm 



(IV-56) 



is the surface area of a critical bubble, or 3/4 it N is the 

^ CO' 
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density of locations (potential sites for bubble nucleation). A potential 
site is any free volume (i.e., volume unoccupied by polymer) in the 
matrix which is large enough to contain a gas molecule. The number of 
such sites can be estimated by assuming that their size follows a 
Gaussian distribution, and that the sura of the volume of all such sites 
equals the total free volume of the material (this assumes that the 
total volume of sites smaller than one gas molecule is negligible). 
Using nitrogen in polystyrene as a typical example, the volume of one 

molecule is 6.59 X lO'^'* ml and the volume of one mer of poly- 
styrene, 1.32 X 10'^^ ml. Since the density of polystyrene is about 
one gm/ml, the total number of sites, N , is approximately 7 X 10^^ sites/ 
cm . This compares favorably with the value of 4 X 10^^ cavities/ml 
for hydrogen in polyethylene found by Frisch (1965). 

Unfortunately, when the program is run with the values typical of 
those which produce foam In the laboratory, the rate of nucleation pre- 
dicted is negligible. For nucleation to be appreciable, the change 
in Gibb's free- energy between the initial. state and the state with a 
cell of critical radius must not be greater than about 60 feT. This was 
shown by the program results and mentioned by Russell (1980). In the 
laboratory, polystyrene was foamed at lU'^C (390''K) after having been 
exposed to 2.4 X 10^ dynes/cm^ (350 psi) gas for 24 hrs. The con- 
centration of gas in the polystyrene is not greater than 0.0306 moles/ cm^ 
(0.00187 moles/in ) which is the saturation concentration at 2.4 X lO'' 
dynes/cm . The dynamic modulus of the material at 393*'C Is 158 psi 
(Waldman, 1980 ) (Tested on the same run of extruded material as was 
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used for the experiments reported in this work). According to experimen- 
tal and theoretical work by Wu (1970), the surface °nergy co- 
efficient for polystyrene is 33 dynes/cm (1.9 X 10"^ lb/in). However, 
looking at the program results, nucleati on is predicted insignificant 
for concentrations equal to or less than 4.3 X 10"^ moles/cm^, even with 
a modulus of 0 and surface tension as low as 7 dynes/cm (one-fifth of that 
of polystyrene). The presence of rubber particles lowers the surface 
energy. Rubber, Itself, has a surface tension of 25 dynes/cm (calculated b> 
the parachor method in Appendix A), so the interfacial tension is 8 
dynes/cm. The energy stored in this interface contributes to the forma- 
tion of cells forming there. In the extreme, if the cell could form 
Its surface entirely from interfacial surface, the value for surface 
tension which should be used in the calculation is 17 dynes/cm 
(8.46 X 10"^ Ib/in) which is still more than twice as high as the 
value predicted acceptable. 

Since the prediction is so far in error, it suggests a flaw in 
the theory. It may be that the classical nurleation theory does not 
apply to this case. Larger radii are associated with slower nucleati on 
rates, but even for rates so slow that nucleatlon is not appreciable 

(the change of Gibb's free energy > 60 feT), the critical radii are 
" -7 

smaller than 25 A (1 X 10 in). A lower bound to the average distance 
between polymer chains was estimated to be about 7 A (see page 65). 
Since the size of the voids in between the polymer molecules follows a 
random distribution. It Is conceivable that voids exist which are the 
size of the critical radius. The actual distribution of cavities 
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as a function of size should be dotermlned to see If this argument Is 
plausible. Since this requires extensive study of the conformational 
statistics of polymer chains, it Is beyond the scope of this work. It 
is also questionable whether the concept of surface energy can still 
be employed, or whether Gibb's free energy must likewise be calculated 
from statistical considerations, 
b) Growth by Diffusion 

Once they are nucleated, cells grow by diffusion from the matrix. 
It was assumed that each cell is surrounded by a diffusion zone across 
which no gas flows, that all diffusion zones are the same size, and that 
they are arranged in a close-packed configuration. The size of the zones 
is determined only by cell density. 

This mode has several limitations which restrict it to short 
diffusion times. First of all, it does not take into account the fact 
th&t larger cells have lower internal pressures and thus higher potential 
for diffusion (Hobbs, 197^. Also, it does not account for the move- 
ment of gas not in a diffusion zone, and it makes whe assumption that 
the cells are regularly spaced. 

The concentration distribution Is found by applying the diffusion 
equation to the diffusion zone. Since the region was spherical and 
uniform, advantage can be taken of the symettry to drop the angular 
terms, leaving 



3x = 
dt 



„ 3^x . 2 3x 
' ;7 * rSF 



(lV-57) 
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The boundary conditions for the zone are that the concentration at the 
edge of the cell is the solubility, K, times the pressure in the cell. 
The restriction that no gas can cross the outside boundary of the zone 
imposes the condition that the concentration gradient be zero at that 
boundary. The initial condition is the concentration gradient shown in 
Figure IV-9. The concentration has the initial value, C^^, everywhere 
except in the depleted region, as shown. 

This was solved with a series solution (Hildebrandt, 1964). The 
particular solutions were of the form: 



and orthogonal in the region of interest. X was found from the character- 
istic equation 




-xV 



X(r,e) = 



Xr 



-tan(XR)cos(Xr) + sin(Xr) 



(IV-58) 



tanXR = 



tan(Xr^) - Xr 



(IV-59) 



Xr,tanlXr,) + 1 



The coefficients. A,, were found from 




R 



(IV-60) 



where 



10 



r2 
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Figure IV-5 



Critical Radius Vs. Concentration (Log Plot) 
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Figure IV-6 Critical Radius Vs. Modulus 
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Figure IV-8 Diagram of Depleted Regions and Diffusion Zones 



Figure IV-9 Concentration Gradient 
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Figure IV-10 Particular Solutions to Diffusion Equation when T - 
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t? - ~tan(XR)cosXr + sIn Xr ,,,, 
•^X " Xr — UV-61i 

For typical values, the first four terms at t = 0 are shown in 
Fifijre IV-IC. 

c) A Computer Model 

A model was developed which combines nucleation rate with diffu- 
sional growth to determine the cell size distribution. It was programmed 
on the MULTICS system at MIT in APL. 

For successive time increments, the program calculates the 
number of cells nucleated and the diffusional growth which occurs during 
each increment. Since the cells first nucleated will grow largest (and 
since the amount of diffusion depends on the radius of the cell), the 
amount of growth Is calculated separately for each time increment, the gas 
concentration is adjusted to account for the nucleation and growth. The 
program and a detailed discussion of it are found in Appendix C. 

The program uses search routines to determine the X's in the dif- 
fusion equation. Unfortunately, after the program was completed In API, 
It was found that the number of X's needed to obtain <:uff1cient1y accurate 
concentration profile? was hig** (at least 12). Since the search routines 
arelengthy {'^ 3 hours), obtaining the cell size distribution for one set 
of conditions would take the batter p»rt of a weekend. Neither the time 
nor the computer funds were available for a complete set of data points 
calculated this way. Although the work could be reprogrammed in a more 
efficient language, enough had been learned about the calculation to 
enable making some additional assumptions and doing part of the calcula* 
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tion by hand. The assumptions and calculation methods will be desribed, 
foTlowed by a description of the results. 

Before delving into that topic, however, it is Important to realize 
that there are three stages during foam formation and growth: nucleation, 
diffusional growth, and matrix relaxation. Nucleation is the transforma- 
tion of small clusters of gas molecules to energetically stable pockets 
of molecules with distinct call walls. The second stage is that of dif- 
fusional growth. During this period the bulk of the gas moves out of the 
matrix and into the cells until equilibrium exists between the concentra- 
tion of gas in the matrix and the internal pressure in the cells. The 
third stage is growth by relaxation of the matrix. Relaxational growth 
allows the material around the cells to lose its orientation, whereas 
nucleation and deformational growth both impose, orientation. 

The program follows foam formation through the first two periods. 
It determines the number of cells nucleated, then calculates the amount 
of gas which diffuses to eacl» cell. It is Irjportant to gain an accurate 
estimate of the total number of cells nucleated to be able to predict 
their final size. 

The program uses classical nucleation theory to predict the cell 
nucleation rate from concentration and other system parameters (-ee 
Section Iv.B.sV That theory shows that the ultimate cell density Is 
limited by the gas concentration in the matrix. Only when the concentra- 
tion exceeds a certain minimum value will nucleation occur. It was as- 
sumed that nucleation does occur, even when the concentration exceeds the 
limit only locally. Nucleated cells deplete the local region of gas and 
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attract gas by diffusion. Between the cells, the concentration will 
drop with time. Thus, the ultimate cell density depends heavily on both 
the Gibb's free energy barrier and the. rate of diffusion of the gas 
through the matrix. 

Three regimes of .foaming behavior can be established by comparing 
the nucleation rate with the rate of diffusion. When the nucleation rate 
is slow (less than 2 X 10^^ cells/inli«),+he few cells which nucleate are 
widespread, and it takes a long time for diffusion to bring the system 
into equilibrium. By the time equilibrium is reached, relaxation of 
the matrix is already significant. Such low nucleation rates produce 
foanis of poor quality with large, widely spaced cells. If the nucleation 
rate is high, cells will develop sofastthat the diffusion length 
( /n where t is the duration of nucleation) is less than the thick- 
ness of the depleted region. In this case, nucleation is halted when 
depleted regions interfere with each other. For intermediate nucleation 
rates (greater than 2 X 10^ and less than 2 X 10^^ cells/in^ sec for 
the polystyrene - Ng system studied) diffusion is of the same order as 
the rate of nucleation. The diffusion process Is well underway by the 
time nucleation stops. In fact, nucleation halts because diffusion 
significantly lowers the gas concentration in the matrix. 

The program is most successful in the middle regime. For high or 
low rates of nucleation, simplifying assumptions make the computer pro- 
gram inapplicable. 

Since the program does not make a detailed analysis of tht conctn- 
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tration distribution in the depleted region, its estimate of the final 
cell density In quickly nucleated foams is dictated by an arbitrary choice 
of concentration distribution. A simpler and more consistent method for 
calculating the ultimate cell densities in these foams is to assume that 
nucleati on halts once the depleted regions interfere with one another. 
This can be calculated easily by hand. 

The total number of cells nucleated can be estimated by assuming 
a close packed configuration of spheres the size of the depleted region. 
Using the argument as proposed in Appendix 6 to calculate the distance 
between closepacked cells. 



where here, r^. Is the radius of the depleted region. For one unit 
volume of completely nucleated foam, the gaseous volume Is 



where R is the cell radius. The remaining volume Is the matrix, so 



p cell - 



(IV-62) 



4/2 r^ 



p cell • |irR^ 



(IV-63) 




p cell • J irR" 



4irR 



,3 



(IV-64) 
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Recalling that the depleted volume is the same as the bubble volume. 



= 2 (IV.65) 



and 



= 1 7^— (IV-66) 

6 / 2 TT 



Rapid nucleation, therefore, will produce 



1 



4/7 r^ 



1 - 



2 ir 



•-d 



per original unit volume of polymer. A typical nucleation radius for the 
fast regime is 5.08 X 10"^ p (2 X 10"® in). In this case, the cell 
density will be 4.15 X 10^^ cells/cm^ sec (6.8 X 10^^ cells/in^ sec). 

For low nucleation rates, the computer model does not take relaxa- 
tion of the plastic into account,. and therefore will over estimate the 
cell density. This matter could be corrected by incorporating a 
modulus which changes with time. For these systems, however, cell 
growth is such a major factor that it would also be necessary to keep 
track of cell "pirating." This is the phenomena in which large cells 
draw gas from smaller ones by diffusion. Since the larger cells have 
lower internal pressures, a concentration gradient is set up between 
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them and the smaller cells. To account for this, it is necessary to 
abandon the assumption that all cells have identically sized diffusion 
zones and allow the size of the diffusion zone to vary with cell size and 
the local configuration of cells (each having a different diameter). This 
introduces a great deal of complication into the program. Making an 
intelligent assumption of a "typical" local cell configuration can 
be tricky. Once It has been established, the size of the diffusion zones 
can be obtained from the condition that the boundary concentrations of 
adjacent regions must match. Hopefully, it would be possible to conceive 
of a packing configuration so that the diffusion regions could be assumed 
spherical, otherwise the calculation would become even more hopelessly 
complicated. 

Since the foams resulting from the slow nucleation regime will have 
large cell sizes and consequently poor quality, a detailed accounting of 
them was not considered worth the effort. Instead, the conditions for 
which this regime occurs have been mapped in Figure IV-11 and IV-12. 

The computer programs is useful for predicting the qualities of 
foam in the intermediate nucleation regime. The program is described in 
full in Appendix C. Because of the extensive calculation time already 
mentioned, the rates of nucleation for various conditions and diffusion 
rates were calculated independently, then were cooruinated by hand. 

The subroutine (or function) which finds the nucleation rate Is 
PROG. In accordance with the classical nucleation theory discussed 
earlier, the change of Gibb's free energy (&G) for bubbles of different 
sizes is determined. The maximum AG is the activation energy which 
must be overcome to produce a stable bubble nucleus. The rate of nuclea- 
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To read Figures IV- 11 to IV~13. 

The nucleation regime is determined by four variables. To de- 
termine if a state is in the fast nucleation regime, pick the graph in 
accordance with the system's temperature, then find the point correspond- 
ing to the surface tension and modulus. The concentration must be above 
that of the equi-concentration line running through that point. Use 
the same procedure for slow nucleation. except the concentration must be 
below that of the equi-concentration line running through that point. 
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Figure IV-12 Bounds Between Fast and Intermediate Diffusion 

Zones - Continued 
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tion follows an Arrehenius relationship. The cell size for which AG 
Is a maximum Is called the critical radius, and It Is assumed that newly 
nucleated cells are this size. The rate of nucleation depends upon the 
size of the diffusion zone, the size of the cell (actually the pressure 
within It), and, of course, the diffusion coefficient. The size of the 
diffusion zone is calculated from the cell density by assuming that the 
zones are all the same size and closepacked. A method to calculate the 
size of the zone is shown in Appendix B. As mentioned before, using dif- 
fusion zones which are all the same size predicts cells of only one size. 
This Is acceptable, however, since in this regime, variations in 
nucleus size are less than 5%. As discussed in Section IV.B.b, the 
diffusion equation was solved for the conditions that the radial derivi- 
tlve of the radial derivitive of the concentration distribution at the 
outside boundary Is zero, and that the concentration at the edge of the 
bubble Is 

C = pressure In bubble ^ jY^ggj 

For the Initial condition, the concentration profile shown In Figure IV-9 
was used. The sloped region (r - R to r^) represents the region depleted 
during nucleation. As an easy and crude approximation. It was assumed to 
be 1/10 of the' outside radius of the depleted zone. Actually, Its thick- 
ness Is .288 X R. 

Once tables of diffusion and nucleation rates were obtained, itera- 
tions to calculate the final cell density were made by hand. An initial 
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time step was chosen so that at least 1% of the final number of cells 
would be produced in subsequent time steps. The later time steps were 
generally somewhat shorter to gain an accurate estimate of the cell 
density. The number of cells produced per unit volume during each time 
step was added until the increase was insignificant. As long as the 
above procedure was followed, variations in time step did not significantly 
alter the results. 

Graphs of cell density versus critical parameters for this nuclea- 
tion reglmeare shown in Figure IV-14, IV-15, IV-16, and IV-17. 

Once an estimate has been made of the number of cells nucleated, 
diffusional and relaxational growth can be predicted. 

Diffusional growth stops once pressure equilibrium is attained 
between the cells and the gas in the surrounding matrix. The size of the 
cells at this point can be calculated by determining the excess of gas 
between the initial concentration and the final equilibrium, then divid- 
ing by the number of cells per unit volume. If the concentration Is In 
moles, this can be converted to cell volume (or radius) by the ideal gas 
law. However, because the pressure of the cell (and hence the equilibrium 
pressure) Is dependent on the cell size, the expression for R Is a quartic 
equation: 

1 = 2 1 f ' 3 (jy^ggj 



This was solved by computer (the program 1s in Appendix D). The resulting 
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Figure lV-15 Cell Density Vs. Modulus 




Figure IV-16 Cell Density Vs. Temperature 
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Figure IV-17 Cell Density Vs. Surface Tension 
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Figure IV-18 Cell Radius VS. Tempbjmure 
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cell sizes are shown in Figures IV- i8 to iv-20. 

The above expression can be used to calculate the relaxational growth 
as well, simply by allowing the modulus to change with time. Graphs of 
cell size versus time (or modulus) are in Figure IV-22 to Figure IV-25. 

This model, as every other, has its limitations. These must be 
recognized and understood before the model can be used with confidence. 

One of the most important limitations has already been discussed 
at length. For the polymer system described, the critical radius pre- 
dicted by classical nucleation theory has sub-atomic dimensions. This 
indicates that a non-classical theory of nucleation must be employed and 
that perhaps even. the concept of nucleation rate is questionable. 

If the program is considered despite this flZiw or applied to systems 
whose critical cell size does not fall within this regime, there are 
other limitations to be realized. One important assumption is that the 
reaction occurs isothermal ly. Because of the speed of both nucleation and 
diffusion, the gas is rapidly decompressed. Since the gas must be re- 
heated, and since the material around the edge of the cell may be cooled 
(thereby experiencing an increase in modulus), the energy required to 
foam the plastic is underestimated by the model. Also, the rapid ex- 
pansion sheds doubt on the expression for Gibb's free energy used in the 
nucleation theory. The expression does not account for momentum effects, 
or the energy needed to overcome increases in modulus due to high strain 
rates. Again, the prediction of the energy needed for nucleation is 
too low. 

In addition to the shortcomings of the nucleation theory, the 
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Figure IV-20 Cell Radius Vs. Surface Tension 
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theory describing the diffusional growth deviates from nature by not ac- 
counting for the change of bubble pressure as the bubble grows. To take 
this, into account would result in a differential equation different from 
that of Fourier, and one more difficult to solve. The approximation used 
gives a lower estimate of the size of the cell, since as the bubble grows* 
the pressure in the cell decreases, allowing more gas to be driven into 
the cell. Also, diffusion would occur at a faster rate than predicted. 

The discrepancies .just described are those between the theoretical 
bases of the model and nature. In addition to these, there are also 
simplifications introduced to be able to express the theory in a computer 
program. 

As already discussed, the program assumes that the cells reorganize 
themselves into a close packed configuration. In addition, the model 
assumes that there are always enough sites to sustain nucleation at the 
classically predicted rate. As with the quasi-steady assumptions, this 
will also predict a value of Gibb's free energy that is too low, and 
thus over estimate the cell density. 
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V. RESULTS 

Temperature, modulus, concentration, solubility, surface energy 

« 

and diffusion coefficient: all of these play major roles in the growth 
Qf cellular foam. To design successful foaming processes, it is important 
to understand the role of each. 

A. Temperature and Modulus 
a) Nucleation 

Since the energy of mixing is proportional to RT In X/Xg, high 
temperatures during nucleation yield a high energy of gas supersaturatlon. 
This leads to high densities of small cells, as shown in Figure IV-15. 
Sin-e the modulus is low when temperature is high, little energy is spent 
deforming the matrix, and the trend of increasing densities of smaller 
cells with higher temperatures is intensified (cell density versus 
modulus for constant temperature Is shown In Figure IV-16, for moduli 
which depend upon temperature, see Figure V-1). It should be noted, how- 
ever, that increasing temperature also causes an increase In the diffu- 
sion coefficient. For the slow and intermediate nucleation rate regimes, 
higher diffusion rates will hasten cell growth and thus lower density. 
For foams produced In the laboratory, all of the material was heated from 
room temperature to above Tg. Thus, the temperature of the matrix is 
increasing during nucleation and growth. Since the gas concentration 
and material characteristics remained unchanged while the oven tempera- 
ture was varied, modulus and surface tension each depended only on 
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Figure V-1 Foam Characteristics for Material with 
Temperature-dependlw Modulus (part a) 
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Figure V-1 Foam Characteristics for Material with 

TEMPERATURE-DEPENiefT fbOUUIS (PART c) 



-109 - 



temperature. Thus, nucleation should have begun at the same critical 
temperature for each sample, regardless of the final temperature. Since 
nucleation is completed in such a short time, the oven temperature is 
expected to influence foam formation most during the relaxational growth 
regime. Inspection of Figure III-7 substantiates this suggestion; the 
size of the cells increases with oven temperature. Since growth was 
halted in these systems at 7.5 minutes, the modulus was not zero, but 
varied with jemperature (see Figure V-3). This would lead to cell sizes 
which increase with higher temperatures. For the intermediate systems 
analyzed by computer, increasing the diffusion coefficient by half gave 
increases in cell density of only 10%. Here the effect of the increase 
of modulus with temperature is much more dramatic. 

b) Diffusional Growth 

During diffusional growth, temperature plays its predominant role 
in changing the diffusion coefficient. High temperature, and thus rapid 
diffusion is responsible for shorter diffusional growth time. In slowly 
nucleated systems which have non-uniform cell sizes, high temperature 
causes larger spread in cell size distribution. In these cases, high 
temperatures are not desirable during diffusional growth since they lead 
to larger maximum cell size. 

c) Kelaxational Growth 

In this realm, high temperatures, by causing the modulus to be 
low, allow bubbles to grow faster and larger. Orientation is destroyed 
more readily. For foams in the fast and intermediate nucleation regimes, 
the effects of temperature and modulus during nucleation outweigh their 
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effect during relaxational growth (see Figure V-2). 

B. Initial Concentration 
a) Nucleation 

Concentration plays its most important role during nucleation. 
With higher supersaturation of gas in the matrix, more energy is avail 
able for forming cells and more cells form. Since the residual super- 
saturition of gas remaining after nucleation is divided among them, a 
higher cell density means smaller cells. Since the amount of residual 
gas has little to do with the initial concentration, realizing the initi 
gas concentration increases the cell density, decreases cell size, and 
improves material properties. 

b; Diffusional and relaxational growth 

The increase in cell size during the diffusional growth is based 
on the excess of gas and the cell density. Since the nucleation rate 
depends quite steeply upon concentration, increasing the initial gas 
concentration decreases cell size, as shown in Figure iv-21. The de- 
pendence of cell size on concentration in laboratory produced foam is 
shown in Figure III-l. 

C. Surface Energy 
a) Nucleation 

Energy must be invested in new surface when cells are formed. 
Thus, the number of cells increases as the surface energy coefficient 
decreases, while their size decreases. See Figure IV-17. 




Figure Cell Density Vs. Surface Tension 
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b) Diffusional Growth 

The surface tension resists growth of the cell and maintains a 
pressure difference between the interior and exterior of the cell. 
Smaller cells have higher internal pressures than larger ones. In the 
slow regime, which has large bubble size variations, this effect causes 
the larger ones to attract gas from the smaller cells. Higher surface 
energy coefficients make this undesirable effect more pronounced, 

c) Relaxational Growth 

Since the surface tension resists cell growth, it aids the mechani 
cal properties during this period. However, as cells grow larger, the 
effects of surface energy become less and less significant when com- 
pared to volumetric quantities, such as modulus. 
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VI. CONCLUSIONS 

Through the course of this work, the following conclusions have 
been reached: 

(1) A method for producing microcellular foam has been demonstrated. 
The material produced has tensile strength which is superior 

to that of conventional foam, and which improves as cell size 
decreases. Hov/ever, the method requires long cycle times 
which eliminate its usefulness as a commercial process. 

(2) Classical nucleation theory was applied to a system consisting 
of a polymer matrix supersaturated with a specified amount of 
gas to determine the rate of nucleation and ultimate cell 
density. The amount of gas diffusing to each cell was esti- 
mated, and finally the growth of the cell due to the relaxation 
of the matrix was determined. Although applying classical 
nucleation theory to polystyrene foams may not be valid, the 
theory predicts that high nucleation rates result from large 
drops in free energy during foaming. Initial and final condi- 
tions which enhance the free energy drop are favorable for the 
production of microcellular foam. These are: 

* Increasing initial gas concentration. 

* Decreasing surface tension. 

* Increasing the temperature of the process. 

* Decreasing the modulus of the materia':. 

(3) Foam has three stages of development: nucleation, diffuslonail 
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growth and growth due to matrix relaxation. Since nucleation 
win halt once the level of gas supersaturation falls below a 
critical value » the type of foam produced depends upon the com- 
petition between nucleation and diffusion for available gas. 
When nucleation occurs so rapidly that gas diffusion is insig- 
nificant, the cell density will be very high and the amount 
of gas which can be absorbed by each cell will be small. 

If diffusion and nucleation have similar time constants, 
diffusion will reduce the gas concentration between the cells 
and ultimately stop nucleation. The foams have intermediate 
cell densities and larger cells,, but the distribution of cell 
sizes is expected to be fairly narrow. 

If nucleation is extremely slow, relaxational growth will 
occur while nucleation and diffusion are taking place. These 
foams should have comparatively poor qualities, with low cell 
densities, large cells, and a wide distribution of cell sizes. 

The number of cells nucleated should be as high as 
possible. To achieve this, nucleation should be completed as 
quickly as possible. If nucleation is triggered during less 
favovable conditions, it does not matter if more favorable 
conditions ensue; large cells have already been nucleated and 
are growing. 

Future Work 

There is still much to learn. In addition to developing a commer- 
cially feasible process for producing microcellular foam, other future 



endeavors are suggested by this work. 

The classical theory of homogeneous nucleation theory does not ap- 
pear to apply to the systems studied. The critical nucleii are predicted 
to be so small that they are comparable to free volume voids existing in 
the polymer. Nucleation should bs examined with a more detailed under- 
sta-'ding of the fluctuations and movements of the polymer chains. Also, 
the event occurs so rapidly that the equilibrium assumptions should be 
questioned. The model should account for the drop in the temperature of 
the material due to the decompression of the gas. 

More data is needed to gain a better quantitative knowledge of the 
effects of various process parameters on the morphology and the strength 
of the foam produced. Waldman (1980) has already begun work on measuring 
the impact strength of the foam. Fracture of foam should be studied, pro- 
bably using the concept of J-integrals. Also, it would be instructive to 
model the amount of orientation imposed during the formation of a bubble, 
how much of it is lost before the foam is cooled, and the effects of such 
orientation on strength. 

Finally, this foaming method should be tried with other combina- 
tions of polymer and gas and adopted for thermo- and chemosets. 

By understanding polymer foams more clearly, they may become im- 
portant engineering materials, offering improved materials properties and 
lower densities. 
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Appendix A 

SURFACE TENSION OF POLYBUTADIENE AT 120«C 

The parachor. P. is a quantity which relates surface tension to 
molecular volume. V (or. in this case, mer volunie) by 



1/4 

(Ap. A.l) 



Values for the parachor can be found by adding contributions fron, each 
atomic group in the .er. A list of the contributions of various groups 
of atoms have been published (Reid 1977). For polybutadiene. the 

parachor is 145 fe- X ^ ] tk 

3 ^ cm^ J * volume of 25'*C 

is 60.87 cm /mole, and the thermal expansion coefficient. 22 X lO'^ 
For a 950c temperature rise, the molecular volume increases to a value of 
64.77 ca /mole. Using the above formula, the surface tension of 120»C is 
25.12 dynes/cm (1.43 X 10"* lb/in). ' 
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Appendix B 
DETERMINATION OF r . (OUTR) 



It was assumed that each bubble is surrounded by a spherical 
"diffusion region" (a control volume which gas neither enters nor 
leaves), and that all the regions are the same size, regardless of the 
size of the cell within. It was also assumed that the regions formed 
a close-packed configuration. For a standard unit cube of a close-packed 
matrix, there is a sphere (radius r^) at each corner and the center of 
each face, so 



Since there are 4 cells per cube (which has iide-length 2a), the cell 
density- is 





(Ap. B.l) 



/Tr^ = 



(Ap. B.2) 



p cell = 




(Ap. B.3) 



Substituting A for a. 
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p cell = 



r \ 

r2 



^ J 



4/2 



(Ap. B.4) 



1 



5.657 r: 



In the program, r^. is called OUTR. 
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Appendix C 

COMPUTER CALCULATION OF CELL SIZE DISTRIBUTION 

This appendix describes the computer program designed to calculate 
the cell size distribution corresponding to specified processing or 
material parameters (modulus, temperature, surface tension, initial gas 
concentration, diffusion coefficient, solubility, etc.). The main pro- 
gram is described first, followed by the subroutines. The subroutine 
scheme is shown in Figure A-C-1, and the program is listed starting 
on page 134. The program was written in APL on the MIT MULTICS system. 

A. Main Program - CELDIST 

1. Input Section 

The first section of the program simply establishes the input 
parameters, and pre-sets certain incrementing values. The conditions are 
printed before the distribution calculation Is made. Ther version shown 
increments through modulus and temperature, calculating the cell size 
distribution for each set of parameters. 

2. Nucleation Section 

The purpose of the second section Is to determine the number of 
cells nucleated during each time Increment. First, however, it determines 
the concentration at which the concentration at which the rate of nucle.i- 
would be negligible (CRITCON: found by the subroutine CUTCONC). so that 
a test can be made to see if nucleation will occur (line 45). SifcomHy. 
It chooses a time step of a fraction the time it would taki? to lUmUU 
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Subroutine Scheme 



Main Program 



CELDIST 



1st Sublevel 2nd Sublevel 3rd Sublevel 
CUTONC 



CRITRAD 



GIBBS 



RATE 



PROG 


CRITRAD 






GIBBS 






RATE 




DIFF 


LAMDA 






6ASIN 






AN 


STRATE 
SLOPE 


NEWSIZE 


















DONELY 


GASIN 






NEUSIZE 





CONCEN 
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the volume of the matrix. (This is chosen because of the gas-volume 
condition described previously.) Finally, the number of cells nucleated 
during that time step is determined. PROG is the function which finds 
the njcleation rate as well as the critical radius, R and the change of 
Gibb's free energy, G. The number of cells nucleated is found by multi- 
plying the rate by the length of the time increment (line 84). A count 
is made (line 75) to be able to keep track of the nucleation groups. 
(This will be used later to index a matrix containing the pertinent 
information about each group.) 

3. Diffusional Growth Section 

The growth of the cells by diffusion is calculated using the 
solution determined in Chapter IV, Section B. 

To determine the growth of the cells by diffusion, each cell was 

assumed to be in a spherical region across which no gas flows. These 

regions are assumed to be close-packed, and their radius is calculated 

by lines 93 and 94. TOTCEL is the total number of cells per cubic thou. 

From the arguments in Appendix B, the radius of the diffusion sphere, 

OUTR, is (1/4 rr TOTCEL)^^^ . The solution found to the diffusion 

equation is the sum of a series of partial solutions, each haying a 

different separation constant X , and summing coefficients A, . The 

" '^n 
details of finding the X's and the A^^'s and incorporating them all into 

a summed solution of the flow of gas (M) into the cell is all included 

in the function (subroutine) DIFF (line 110), which depends on the number 

of the group and the size of the cell. The codes between lines 93 and 110 

establish the concentration gradient surrounding newly nucleated cells. 
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This gradient is shown in Figure IV-9 it ic tHui^^A m 

3UIC iv ^. It IS divided into 10 segments; 

those m the depleted region are attributed a zero CONC (line 105) 
(PRESS/K, the actual concentration at the boundary of the cell will be 
added to. the final solution.) The segment containing the boundary of 
the depleted zone is slanted between 0 and the initial concentration (so 
its slope is initial CONC/seg«ent length. ieg«„t length is callea CHRI 
In the pr-^ram (see line 96). The remaining segments are at the Initial 
concentration. This scheme Is facilitated by Initially assigning the 
Initial concentration to all of the points in CONSSAVC (the values in 
CONSAVE represent the concentrations at the end points of each segment 
- this would make 11. The twelfth value is simply used as a place to 
store the segment length. CHRI. since it will be needed for later 
reference. CONSAVE [12] Is reassigned in line 101. The loop (LOUP) 
assigns 0 to points in CONSAVE from the cell radius out until the 
numbers of segments exceeds that which fit into the depleted region 
(line 104 and 98). 

Once DIFF has established the flux of gas into the cell, NEWSIZE 
determines the cell size after all this gas has been added. The first 
part of the section (lines 91 to 117) calculates the diffuslonal growth 
of newly nucleated cells, and lines 120 to 132 calculate the growth 
of previously nucleated cells. For diffusion Into groups of cells which 
are not newly nucleated. DIFF Is also used, but this time, the Initial 
concentration Is that established during the previous time Increment. 
The LOOP which evaluates growth of these cells extends from lines 120 
to 132. The concentration used for calculating the next nucleatlon rate 
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is that at the outside edge of the diffusion boundary of the first group 
of cells. When these calculations are completed a new time is assigned 
(line 121) and the concentration is re-evaluated. Once nucleation" is 
complete, cell growth is determined by DONELY. 

B. Subroutines 

1. PROG 

PROG unites the three functions: CRITRAD, GIBBS and RATE, to 
establish the rate of nucleation, critical radius, and this change of 
Gibb's free energy during nucleation. If CRITRAD does not find a 
critical radius (R = 0), then the computer returns to the main program 
(where it is aborted again). 

2. CUTCONC 

CUTCONC also unites CRITRAD, GIBBS and RATE, but this time Incre- 
ments DOCON from the Initial concentration to determine the concentration 
at which nucleation will become negligible 

3. DIFF 

DIFF calculates the flow of gas into the cell (in moles per cell) 
by calling three other functions, LAMDA, EQCOEF and GASIN. As already 
mentioned, the flux of gas into the cell is calculated by solving diffu- 
sion equation for a spherical region surrounding the cell with Insulated 
boundary conditions, then applying Fourier's equation at the cell 
boundary. The basis for the Insular condition Is the assumption that 
gas will diffuse only to the closest cell, and therefore, there are 
boundaries between the cells across which flow will not occur. The pro- 
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this IS not true for cells nf w4rf.i. • 

ceils of widely varying radius, where the boundary 

.s «c Closer to the s.a„er cell. :f .he discrepancy U Ur,e enough. 

s will actually leave the seller cell and .0 to the larger one 
This ,s discussed In Chapter V. The subroutine calculates the first 8 
value. Of X (line 1). .en the values of . (1,. , and ,0, dete^lnes 

he f ux Of gas Into the cell (line 1., and the new cell si.e (line ,4,. 

■Z " ' ^WT. and the re- 

sulting concentration by calling SAVECONC. 

4. NEWSIZE 

NEWSIZE uses the Ideal gas law tp approximate the »,ount of gas 

flowed into the cell during the last f,-™. • 

, . ^ increment. Finally, the Ideal 

gas law u used again to detemine the final size of the cell. 
5. CONCEN 

CONCEN dete^lnes the concentration at any radius In the diffusant 
-gion. It does this by getting the values of X and A fro™ EQMAT for 

e particular cell gr^p ,„est1on. and plugging the™ Into the solu- 
tion for concentration. Again, the solution Is a series solution, so the 
progrm loops until four tenns have been su^ed. 
6. DONELV 

DONELV calculates the growth of the cells after nuCeatlon has 
ceased. It takes the values of X and A established during the previous 
tl^e increment and detennlnes the growth of each group of cells as t1^ 
P-Sresses by taking gas In. This procedure leaves the c«,centrat.o„ 
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distributions intact from one time increment to the next. 
7. CRITRAD 

CRITRAD searches for the radius at which the change of Gibb's free 
energy is maximum, called the critical radius. First, it finds the 
order of magnitude by rising an exponential search, then it changes to 
a decimal search. To find the order of magnitude, it creates the vector 
TRYR. Consisting of a list of exponentially increasing numbers (line 2). 
It then determines the derivitive of the Gibb's free energy change with 
respect to R at each of the radii listed in TRYR. Line 11 establishes 
a vector (HOLD) the same lengch as TRYR, but having a 1 ir the derivitive 
Is negative, and a 0 If it is positive (line 11). Line 12 takes the 
successive sum of the components in HOLD. (O^O^O^l^l^l) -►(0 0 0 1 2 i). 
IND is assigned a vector having a 1 only where the vector has a 1 (and 

0 elsewhere), line 13. By multiplying IND with TRYR. the radius R Is 
found which Is an over approximation to the critical radius. Now, TRYR 
Is reassigned a new vector, this time a list of numbers increasing 
decimally up to the value R. For the exponential search, this Is done 
by assigning log R to START. Then multiplying a list of numbers from 

1 to 10 by 10 raised to (START-1). line 17. (The first term in 17 Is 

0 because court is 1, and 10 raised to START is R. See Gilman and Rose 
(1976) for a PL calculation SYNTAX; it calculates from right to left 
without regard for the type of operation). NUMB is Increased by 1 and the 
calculation is repeated starting from LUPE on line 4. Now, when line 14 
is calculated, R is multidigital. Its exponent Is stored in START (by 
rounding up the logarithm), and now the first term In line 17 decreases 
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the last digit in R by l, then add^ tn 

. ^ * *° «''»'ay of numbers from I to 

'0 in the next smallest diait Thi* 4e 

aign. This IS repeated 8 times for accuracy. 

8. GIBBS 

GIB8S Sl«„,y ccuutes the change ,f G.^b-s free energy once the 

or «c., u .etennlned. There are three tenns 03. D4 and OS; the 

are added together. 
9. RATE 

RATE ,s a straightforward calculation of the rate at which cells 

fona. The rate equation was discussed earlier a„H it w - 

^ , earner, and it was shown that the 

rate of nucleatlon, 

>*er= ,^ is the frequency at which gas «„,ecules strike a surface. S is 
the surface area (calculated ,n lines 1 ,„d 2). ^he n„her of ' 

msm. nucleation sites (6.55 X 10^ sites/.n', and Z is the Zeldonich 
actor (calculated as line 11). u the change in Gibh-s free energy 
for . cell Of critic, si.e and was calculated by GIBBS. K ,s Bolt^nan's 
constant, and T the temperature. 
10. LAMDA 

LAMDA determines the separation coefficients of the diffusion equa- 
'''^ » '"""'te number of these constants, only 
IS are determined by this calculation. The order of „gnitude ,f 
the first X is established by letting TRYL equal a list of 1000 numbers 
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between 0 and 1/Jog OUTR. This choice of TRYL insures that two roots 
do not appear in the list, that the root does not occur between 0 and the 
first value in the list (1/1000 . log OUTR): and that the root is de- 
termined with sufficient accuracy. The characteristic equation that X 
must satisfy is evaluated by line 15. ONEL must be 1 if the equation is 
satisfied. To search for A. ONEL is compared with a vector RONEL (which 
Is ONEL rotated by one element). All intervals over which ONEL changes 
from less than one to greater than one are tested for slope to see if they 
are zeros or poles. The zeros. are recorded as values for X. The calcula- 
tion is repeated for higher values of TRYL (the second set of values from 
l/(log OUTR) to 2/(log OUTR) until 15 X's are obtained. 

Due to the extensive search routines and the time inefficiencies 
of APL, this sub-routine often takes three hours to run on the MIT 
NULTICS system. 

11. NEWEQCO 

Only the first part of NEWEQCO, which detennines the denonenator 
of A, Is used. The numerator Is determined by AN (which must be run 
after NEWEQCO and writes over many of the values found by the later part 
of NEWEQCO). It is a fairly straight forward - albeit long - expression 
of the integration stated for A in Chapter IV, Section B. Lines 7 and 8 
are concerned with the possibility that the diffusion zone may be 
smaller in this time Increment than the last. In this event, the CONC 
distribution is cut off cn the outside. 

12. GASIN 

GASIN plugs the values for X and A into the expression for the 
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derivitive of the concentration at the cell boundary. Multiplying this 
by the diffusion coefficient and by AnR^ determines the amount of gas 
which flows into the cell. 

13. AN 

AN is another long computation, but this program calls on sub- 
routines SLOPE and STRATE to keep things simpler. The program sums the 

contribution of the rectangular and triangular parts of the CONC distribu- 
tion. 

14. SAVECONC 

This program simply loops to recalculate the concentrations stored 
in CONSAVE. {This could be done with a matrix approach as opposed to 
a loop.) The function determining the new concentration is CONCEN 
(line 6). 
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f CELLIST 

MaD^250 

COMJCOf-CO 

SUR ♦ MODf.MOri+250 

•4(MOD>2000>/0 
^^Q4.0» 00301 

COMC4-CO 
TEMP4-390 
MOtiU J TEMP-^TEMF +IO 
COWCf-CO 

^(TEMI=^>420)/SUR: 

nPRXNT OUT IWPUT PARAMETERS 



CELL SIZE BISTRIBUTIOM SatEyUAT J OW ' 



INPUT parameters; " 

I 

COMCEMTRATIONs ' 
CONC 

■MODULUS, TEMPERATURE, SURFACE TEHSZOM, AMD l/SOLUBILXTY= • 
MOD 
TEMP 
3URFT 
K 



At******A****A**AA*« aMUCLEATION SECTION A *AA*AAAAft**AA*AAAAA 

n 

•CALCULATED VAL UES $ • 



hFIRST, find the CONC AT WHICH NUCLEATION WILL 
n NOT OCCUR « 

"THE CONCENTRATION AT WHICH NUCLEATION WILL NOT OCCUR IS • 

CRITCONf-CUTCONC CONC 

CRZTCON 

ft FIND THE INITAL RATE OF NUCLEATION 
PROG 

<G^(1 ♦2218E'"22xTEMP) ) >85)/c^ON 

ANOW, CHOOSE A TIME STEP ^ LET XT FEDEPENDENT ON THC RATB 

SO THAT ONE TENTH OF THE CELLS ARE NUCLEATED DUAlXNa 
^ EACH STEP^ 
•THE TIME STEP IS « 
->(R=0)/CON 

T1ME4-(1^(4,188X(RXJ) ) )X(0*0001X1E-9^<RXR) > 
C TIME > 1 > /TSTOP 
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11543 BACKJTIME 
C553 RUMTIMEf.TI.ME: 

11561 flFIWAUL-V, CAUCULATE THE CELLS NUCLEATED nUPrlNG THAT TIME STEP 

C57D courn>o 

L581 TOTCEU^O 

C593 MATRMf.2 4 rO 

1621 START NUCLEATIOM i 

11633 • • 

C643 EOMAT4-2 2 8 fO 

C653 AGIAWJ • « 

C6^3 • • 
C673 'time:« 

C683 RUNTIME 
C693 'COMCr=« 

C703 coHC 
C713 «=p:og 

C723 'EGUIL. COMCrr' 
C733 EQCOWCf.PRESSC13-rK 

C743 EOCONC 

C753 COUMT^-COUMT^'l 

C763 • ' 

C773 •THE CRXTXCAU RADIUS IS" 

C783 

C793 'THE CHAI-IGE OF GIBBS FREE ENERGY (IH IM-.LB) IS • 

C803 G 

C813 •THE RATE OF MUCLEATION XII CEL.LS/CU ^ THOU/SEC IS* 

CB23 ^ 

C833 'THE MUMBER OF CELLS NUCLEATED FER CU « THOU IS* 

C843 NOi-wlxTlME 

IB51 wo 

C863 PRESS4.PRESSCI3 

C873 A 

C883 iM(iinnttniLfiHkiit*fLfLkAHAHiif[fit^xFFUsiON section* a aaai^ a a a 

C893 A 

C903 • • 

C913 ^CALCULATE THE DIFFUSIONAL GROWTH^ 

C923 TOTCEL<-TOTCEL+NO 

C933 0UTR4- ( ( < 2* ♦ 5 ) ^ ' OTCEL ) * ( l-r3 ) ) X IE-3 

C943 -♦(COUNTp(!l ) /CHEAT 

C953 JCEPR^Rx 1*259921 

1:963 CHRli.(OUTR-R)4^10 
C973 ^(OUTR<R) /OUTDONE 

C983 i:iEFN4.r<»EP'P:-^CHRl) 
C993 CONCr4-CONC-(FRESS4-K ) 

C1003 consaveccount; 3*.concy 

II 1013 CONSAVECCOUNT J 123 *-CHRl 

1:1023 iNr.x<-o 

C1033 LOUFJ INDXf-INDX+1 

C1043 INDK>DEFN)/GOON 

C1053 CONSAVECCOUNT} lND>:3<-0 
riOAT ^LOUP 
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1.1073 GOOHJ •COHSAVErr* 

tllOJ »->xi-K R 

UlllJ -»(M<0)/MODU 

1.1123 PROCEEDi: 

C1133 S1SE4.TIME NEWSISE R 

L1153 veCTOft*.ROMTlME,NO,S12E,R 
tll«3 M«TRX<.MATRX,C13VECTOR 
Lll/J -KCOUMT=l)/NEWTIME 

LllUJ n 

Lliy3 OROOP«-0 

I i^oj MEwe: GROUP- ♦•GROUP+l 

tl2l3 ^( GROUP > COUNT) /NEWT IME 

C1223 

C1233 ho<-matrxcb; 23 

C1243 R«-matrxcb»33 

111253 ^((TIME-(OOTR*2>><lE4)/MO»»Pl 

C1263 M<-GROUP DIFF R 

i:i273 -♦(M<0)/>«ODU 
Cl''83 P'ROCEEi>2: 

C1293 SIZE«-TIME HEi-SlSF R 

L1303 matrxcbJ33*-si2E 

ll^ll outi>one:4(couht«i)/»over 

ri323 -»hewg 

C1333 donuy:i>one;.y eomat 

".353 neS^^mejruhtime^-ri.ntime+time 

l13A3 matrk 

ri373 CONC«.CONSAVECini3+EC»CONC 

C1383 'NEW concentrations* 
CI 39 3 CONC 

ri403 REL.CRTCaN«.CRlTCON+EOCONC 
C1413 H(CONC<CRlTCON)/I>ONLY 
C1423 •♦AGIAN 
C1433 TSTORJSTPM 
C1443 TIME«.,1 

C1453 •»»«'^'' 
C1463 modifi:m«.o 
C1473 -»»^** CEEDl 
C1483 NOl>iF2tM«-0 

[III] cHEATtCONSAVECCOUNT-, 3-CONSAVECl 
vuv,.«;'.X«e STEP TOO LARGE.' 



vr ■ 
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V SAVECOWC GROUP' 
RNO^O 

CHRf. ( OUTR-R ) 4- 1 0 
CHAMGR4-F: 
EMCORIH: { RN0f.RH04*l 
->(P:WO> 11 )/GO 

C0MSAVEi:GR0UP$RN03<-CHANGR CONCEM GROUP 
CHAMGR^CHANGR4-CHR 
^EMCOPrE 
GO i CONSAVEi: GROUP' ; 12Df-CHR 
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N7 R^-COWC CRITRAri PAR AM 

Hr:r"'' '^-^ 

?it'2'7182812«<i>Er.coNc*cOHC>>xTE«P,73.4i,coNCx(TRTR,1)xio SAA/ 
C73 »C0HCt--2xsuRFTx( (TRre,2)xli) '<••<". «l ■<«i>xl^. j646 

C83 02.-DC0HCX < COHC+DEPCOMC ) xTE«Px73 . 41 X ( TRTR,3 , ,4 , , 888 

till HOLD<.tiGIBBS<0 

C123 HOLDf.+\HOUD 

CI 33 lW&<-HOLr.=:l 
C143 IMBxTRYR) 

C153 ->(R=0)/0 

C163 STARTf-r ( 10»f< ) 
C193 YES<-MUMB<10 

C202 -»i.ufexVes 



V 64.CONC GIBBS R 

C 1 3 COHCCEPL^ (14.7+(2XSUP:FT^R) + (M0»J.3,14159))AK 

C23 TOTGAS*.(R»3)xCOMCx4»1888 

C33 »3« <(R:*3)X< 61. 757+1, 333XMOD)) 

C43 ''4<-(R*2)xsuf::FTxl2.566A 

C51 *='5<-<e<COWCDEPL^COWC) )xTEMPx73,61XTOTOAS 

C63 G«.B3+D4+B5 
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V J+.COWC KATE B: 

C103 

C123 -»(G>iooxl,2218E-22xTEi;T;J° 22))*. 5 
C133 ^< (G^d .2218E~22xTEMP) ) >85)/0 
C14D '^3«-*(-G^<l,2218E-22xTEMP) ) 
CI J J -'«-23x2ix£.555E22 
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CID Rf-COWC CRITFtAEi PARAM 

C23 ^(P:=0)/0 

C33 Gf-CONC GIB&S R 

C4D ( ( ( 1 ♦ 2218e:-"22xtemf- ) ) > 85 ) /wew 

C53 RATE R 

CA3 V 



V M4-GROUP DIFF R 

112 LAHB^OUTR LAMDA R 

C23 AVEC^-O 

C33 LAMNOf-O 

lAl STG^^o 

C5] SUMJL-AMWOf.L.AMMO+1 

CA3 ^(UAMNO>8)/GO 

C7D 4( ( (La2)xi)Iffcoxtime) >80)/go 

[83 Lf-L-AMBCLAMNO] 

C9] A^-MEIMEQCO 

ClOD 

cm A 

C123 AVECf.AVEC,A 

C133 M^-R GASZM TIME 

[;i43 SIG4-M4-SIG 

C153 •♦SUM 

CI 6] GO J AVEC<-14.AVEC 

1:17] AAUGH*.UAMB,C ♦Sa^VEC 

C1B3 EOMATf-EOMAT, CIHAAUGH 

L193 SAVECOWC GROUP 

C203 M4.SIG 
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CIJ OUTR: 

1153 L.AMBf.Q 
C63 HONEJV2<.TRrUxOUTR 



C203 L.OUT; .=LAMu^, 

C223 LAM»<.ig^L«„p 
L233 L.AMB 



PGASI».iCQ3^ 
V Mf-R GASIM T 

X2a c:i<-2o(LxR> 

rl^ I^''''*-^l-<(30(i-xR))xci) 

•< '-''<'^*2)) > + (TeftM2T.R) 



AW 



1:33 

C43 



1:53 ^ «-<s._AH;^3;:Ax?,:i:;:;>;<^'^'^«^^ (chr^r)), 



V Af.MEWEOCO 

^ SUMMf-O 

3 t(EWSLOF>E«-0 

3 TO3T<.0 

3 ^AVEMUMf-0 

3 RWOt-O 

3 OUTP: 

3 HISS<.L< <OUTR-R)^COW5C12J) 

3 APQUTR«.(HISSxCOMSC12Il 

3 'AF-PROXIMATION OF OUTR= ' 

|0D AP-OUTR 

J13 si«.io<LxR) 

422 S2<-l0(t-XAF-0UTR) 

133 cif.2o(LxR) 

'H3 C2<-20<LXAP0UTR) 

ISl T1<-30(UXR) 

J^3 TERM1«.< AP'OUTR-R)~2 

173 TERM24-< <C2XS2)-(C1XS1))-S.(2XU> 

IQD TEU-(T1*2)X<TERM1+TERM2) 

193 T2<-(TERM1-TERM2) 

203 T3<.Tix<<S2*2)-(Sl*2)>-5-<L.x2) 

E213 T3<-T3x2 

f223 fENOM<-(TEl + T2-T3)-r(L*2) 

£233 'r<ENOM=' 

C243 DEMOM 

[253 cHR<.coMSC123 

C263 CHAN6R4.R 

•[273 lNTeGJRNO«.ftNO+l 

i[283 -»(RNO>HISS)/OM 

C293 -♦(RMO> ID/ON 

C303 CRl«.(CHRx(RNO-l) )+R 

C313 CR2<-<CHRxRNO)+R 

[323 cc<.2G(i-xCRi) 

C33J sr*.ia<i-xCRi) 

[343 CR<.2o(«-xCR2) 

[353 SR<.ici<LxCR2) 

[363 supiA*.<-( (CRi^i.)*2)xcc) + ( < (2xcri >^<l*3 ) > xsc) + ( ( (2-f<L.*4) ) xcc) ) 

[373 SLP2A<-< < (CRir-L )*2)xsc) + < ( <2xcri)^(L*3) )xcc>-( <2T-(i.ft4) )xsc) 

[383 si.F iB<.<-( (CR2-ri.)*2)xCR) + < < ( 2x CR2 > -r < I- *3 ) ) X SR ) + ( ( < 2-^- < »- *4 ) > XCR ) ) 

C393 Sl.F-2B<-( < (CR2'«-L)*2)XSR) + ( ((2XCR2)-f<L.*3) )XCR)-( (2+<L.«4> >XSR) 

[403 SLOI=-ERTi4-(lTCHR)x(SI-FlA-< (T1XSL.P2A) ) > 

[413 SLOPEPT24-<l-5-CHR)X<SlLRlB-((TlxSUP2f ) ) ) 

[423 SLOPEI^T«.(SUOPEPT2-SLOPEPT1)X (COMSCRWO+13-COMSCRW03) 

[433 STR1«.< (CR2*l.)XCR)-( ( (l-rL)*2>XSR) 

[443 STR2+-( <CR2-i.l-)xSR) + (< <1^L)*2>XCR) 

[453 STR3*.< (CRl^L)xCC)-.( ( <1^I_)*2)X(10(I-X( (CHRx<RWO-l ) )^R) > > ) 

[463 STR4«.< <CRi^i.)xsc) + ( ( <l4.i.)*2>xcc) 

[473 STRAIGHT«.<l^l.)X( (-(STR1+T1XSTR2> ) + (STR3+Tlx6TR4> ) 

[483 ^<RMO<2>/EGG 

[493 TOST«.T ST+STRAIGHT 

11503 EGG J 

[SI 3 .•«UMv-(Sl.0PEPT+C0MSCRM03xSTRAIGHT )+E«KHOM 

[523 SUMMf-SUMM-f NUM 

C533 CHAMGRf.CHANGR'l-CHR 
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V STRAIGHT<.STf?:ATE RVAU 

C4D STR2f. ( ( RVAU^L ) XSR ) , < ( ( 1-L. ) *2 ) XCR ; 
CoD l^L) X (STRl+Tl XSTRO) 



V SLO»>EPTf.Sl.OPE RVAL. 
C13 SRf.lo(j_xRVAL) 

C23 C-«-2o<i-xRVAL> 
C33 ' 

C43 

C53 ^ slopef;;supi^7tIJs^p2) 
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III 3 Mf.MxT 

C23 MOf.(PP:ESSx4»1888x<H:A3) )-^(73*61xTEMF 

i: 4 3 HEW VOLf- < MOLEG X 73 • 6 1 X TEMF ) ~F P:ESS 

C53 S12E*. ( NEWVOL4.4 ♦ 1888 ) * . 3333333 

V 



VDONELtXD] 
V 3&OMELY EGMAT 
HID P:UNSTEFf.l^( (LA2>Xr.IFFCO) 
£21 KOUMTf-O 

C33 REPEAT «i%UNTZMEf.R:UHTIME+RUNSTEP 
1:4] MATRX 

C53 

£61 'TIME}' 

£72 KOUMT^-KOUHT^l 

C83 RUNTIME 

£91 ->(KOUNT>85)/0 

C103 ->( (PcUMTIMEx(L*2)XDlFFCO) >85)/0 

CI 13 GROUPf-O 

CI 23 ftOWL:GROUPf.GROUP+l 

CI 33 GROUP >COUMT) /REPEAT 

C143 LAMNOf-O 

tl53 Bf-GROUP+2 

C163 R*-MATR;;cBf 33 
C173 M<-o 

C183 UAM{L.AMM0*.LAMM0^1 
i:l93 -KLAMMO>8>/GRR 

C203 Lf-EQMATCBf 1 JLAMMO] 
C21 3 A^-EOMATCB J 2 f LAMMO3 

C223 (RUMTiMEx(L*2)x»iPFCo> >a5)/0 

C233 SAGf.R GASXN RUNTIME 

C243 M<.SAG+M 

C253 ->L«M 

C263 GRRJ->(M<0)/POWU 

C273 MATRMCBF33<-f^UMSTEP WEWSIZE R 

C283 ->ROwu 
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V HEWCON^CHOSER COMCEW GROUP 

C2J MEWCOW<.0 
C33 SOLJADf.Af+1 

C53 »*-GROUR+2 

C71 ^<^ECtMATCBJ2f 

C83 MEWC^-TIME C CHOSER 

C93 ^^EWCOW^WEWCOWVWEWC 

C103 -♦SOL; 
V 
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Appendix D 

CELL SIZE AS A FUNCTION OF MODULUS 
(FINAL) 

This program calculates cell size as a function of modulus (or 
as a function of time as the material relates). 



7 FIMAI- 
C2II PETE JMODt-MOD- 100 

C33 -»(MOi><0)/0 

1:43 TRR«-5E-8+<TRx\1000) 

C53 UPP>ER<.73.61xTEMP^(4.1888XTOTCELX(TRR*3) ) 

C6D LOWER<.(COMC^ ( ( SURFT x2-^TRR ) + ( MOD*3 . 14159 ) + 14 . 7 ) )- < 1-5-187000 > 

C73 SC«N«.LOWERxUP'P'ERXlE-9 

CBl BR«.SCAN<1 

C9D DO<.BR\l 

dOD MOD,TRRC»0] 

C113 -»PETE 
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Appendix E 
DIFFUSION AND SOLUBILITY COEFFICIENTS 
FOR Ng GAS IN POLYSTYRENE 

chamber with 2.4 X lo' dynes/cm^ „ „ , 

, ^ "2 9« for varying lengths of time. A 

"eaur analytic balance (Model „ ,s , was .,sed. The change of 
weight d„e to absorption of the gas was measurable with at least two 

J-.ca„t digits (the smallest change was 0.00.8 gms for a piece which 
onginally weigh ed 9.6529 gms [after 6 hours N, at 2.4 X lo' 
^esults Of this are Shown in ngure E.I. Both the solubility of 

he « ,n polystyrene and the diffusion coefficient can be calculated 
Trom this data. 

A. Solubnity 

Assuming Henry's law holds. 

0.0016 gms N, 

K = • . 1 

gijl PS ~ 1^ 

2.4 X ID' dynes/cn.2^ (Ap. E.l) 

= 6..;x,0-" S^/dynes/cm^ 

B. Diffusion . . 

For a Sheet of polymer suddenly i«ersed In a penetrant vapor at a 
Higher pressure, the amount of penetrant entering the sheet after short 
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aoo2o - 



0.0015 - 



^ aoolo 

I 

OjOOOS 




12 24 36 48 
Time in pressure chamber (hrs) 



Figure Ap.E.1 Gas Concentration Vs. Time 
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time periods is M(t). where 



il!(t) _ 4 rot 

W) < IT (Dudek. 1979) ^ 

(Ap. E.2) 



M(^) is the amount of penetrant which enters the sheet at t = « p i 
the diffusion coefficient, and a suitably short ti.e step Is onl for 



which 



/pt ^ \ 

'■ * (Ap. E.3) 



This expressfon fs based on the as- 
sumption that the Sheet was at e,u„1br1« before l^erslon. that the 
diffusion coefficient Is not dependent upon concentration, and that 
equilibrium is reached Instantly at the surfaces of the sheet 

The diffusion coefficient is thus calculated fro. the daU ,„ 
Figure A-E-1 to be 3.C7 Kio-'^crnVjic. 
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Appendix F 

PRODUCTION OF MICROCELLULAR FOAM BY EXTRUSION 



The number of cells nucleated should be as high as possible. To 
thieve this, nucleatUn should be completed as quickly as possible'.* If 
ro'cleation is triggered during less favorable conditions, it does not 
liitter if more favorable conditions ensue; large cells have already been 
uclea ted and are growing. 

F.N. Cogswell (1972) has evaluated the pressure drop which occurs 
across a die for making extruded rod. He assumed that the drop from the 
reservoir behind the die to the atmosphere could be evaluated by. summing 
■three component's: the pressure reduction due to shear in the diel, the 
^rop due to extensional flow through the die, and the drop incurred by 
"the material exiting the reservoir. If a die is assumed which ha% 
la.nds of zero length, then the first two components are zero, and the 
pressure can be calculated as the raaterikl exits th6 reserv6i¥r'%ii^ 
,a function of the distance from the die opening, xo find ' 
out whether or not microcellular foam can be produced by extrusion. It 
is important to know how quickly the pressure drops. In Append4xllll 
of his paper. Cogswell sums over an infinite series of dflntril volumes 
to find the net drop from T = » to T = 0. The ratio of the ^o^zcm^, 
one differential element to the total drop can be evaluated as 
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Where P is pressure, p is the f-n+,i 

IS the total pressure drop, ^ = , ^ , 

denoted as n), and a, = L dn ]3 melt, index (usually 

^-.ve. no« in U a .unction o. t, " "^^"^ " 

-OP . a CO,..,,.,. ^-^'2 " ~ 

to be the shape of the region 



a 



0. 



2 



(Ap. F.2) 



and 



^'here u the radius of the die. 



Where t and u are apparent viscosities (usual 
spectively). 



(Ap. F.3) 
ly termed A and n, re- 



If the flow were Newtonian (which 



shear rates below lO'^/sec), Eq. 



occurs for polystyrene melts at 



Eq. CAp. F. 1) reduces to 



^ « 1 



Substituting for o). 



df> 



4 



t + 



+ dT) + n 



(Ap. F.4) 
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-3 

If the pressure drop is to occur in less than 10 sec, and if the 
material is moving toward the die at 0.05 in/sec, then this equation must 
predict that the pressure drop be negligible 5 X 10"^ in before the die. 
Using an of 0.025 in and setting dH = 1 X 10"^ in (using a dH of 
5 X 10'^ in only affects the 3rd significant digit), integration (done 
on a Tl 59 programmable calculator) shows that the pressure has already 
dropped 99.5?;. Thus, pre-nucleation has occurred. The distance over 
which the drop is significant is of the same order as the die radius. 
Reducing the die opening to 5 X 10"^ is unfeasible. To make the orifice 
a useful size, the melt velocity must be higher. For instance, an orifice 
with a 0.01 in radius requires a melt velocity of 10 in/sec. This is con- 
ceivable, but larger die openings demand prohibitive extrudate speeds. 

If the fluid is not assumed to be Newtonian, examination of 
Eq. (Ap. F.l) shows that the only way to make |^ decrease rapidly with 
increasing T is to make x increase rapidly with T. This depends upon 
the behavior of the function 

1 T 

♦ ■ 771 'rrj 

T 

Figure (Ap. F.2) shows that i* should be small as possible. But 

so * should be as large as possible. For most viscoelastic polymers. 
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However , U Uss than Thus, the Ne„to„„„ nu,d already .„a„zed 
- the best case. That analysis showed that it would be difficult to 
extrude ,„ but very narrow profiles of n,1croce,l„,ar foa™. The only 
hope for larger shapes Is to find a po1yn,er-g.s c«„b1nat,„n which has 
a very ^o. diffusion coefficient, even at extrusion te„peratu,.s. In 
this way. the ll™it of lo'^ sec might be extended. 
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Appendix G 
INTEGRATION TO FIND VISCOUS ENERGY TERM 



^>IIJTEDF:ATECDDV 
9 IKTEGRCiTE 

Z^-} CUErE<-l+( (Ir-R) *3) 

1:43 A<-(CU£<e* (-2-!-3) )-l 

USD CUBE ^ (-5^3) >-^P: 

i:^,] TERMi<.AvE'):8x3. 14159 
1*73 c*. ( CUBE *< 1-3) )-i 

i:&3 E'«-(CUEiE* (-2-!-3) )-rR 

1:93 TERM2«-Cxt.x8x3. 14159 
i:i03 c.um^(+/termi) + (+/tef::M2) 

1:113 SUM<.SUMX .001 



